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ABSTRACT 

GALEV evolutionary synthesis models describe the evolution of stellar populations 
in general, of star clusters as well as of galaxies, both in terms of resolved stellar 
populations and of integrated light properties over cosmological timescales of ^ 13 
Gyr from the onset of star formation shortly after the Big Bang until today. 

For galaxies, GALEV includes a simultaneous treatment of the chemical evolution 
of the gas and the spectral evolution of the stellar content, allowing for what we call 
a chemically consistent treatment: We use input physics (stellar evolutionary tracks, 
stellar yields and model atmospheres) for a large range of metallicities and consistently 
account for the increasing initial abundances of successive stellar generations. 

Here we present the latest version of the GALEV evolutionary synthesis models that 
are now interactively available at http:/ /www. galev.org. We review the currently 
used input physics, and also give details on how this physics is implemented in practice. 
We explain how to use the interactive web- interface to generate models for user-defined 
parameters and also give a range of applications that can be studied using GALEV, 
ranging from star clusters, undisturbed galaxies of various types E . . . Sd to starburst- 
and dwarf galaxies, both in the local and the high-rcdshift universe. 

Key words: stars: evolution - galaxies: evolution - galaxies: formation - galaxies: 
stellar content - galaxies: abundances - galaxies: star clusters 
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1 INTRODUCTION 

GALEV (short for GALaxy EVolution) evolutionary synthe- 
sis models have been developed over many years. They were 
published in several steps and under a variety of first au- 
thor names, reflecting the number of students who have con- 
tributed their respective shares to the development. GALEV 
models include the spectral evolutionary synthesis of a stel- 
lar population with arbitrary star formation history on the 
basis of the time evolution of the stellar population in the 
Hertzsprung-Russell diagram, as well as a detailed chemi- 
cal evolution model for the ISM in terms of a large num- 
ber of individual element abundances. GALEV models have 
a wide range of application from star clusters (SCs) to re- 
solved nearby galaxies, to more distant galaxies observed 
in terms of integrated spectra and photometry, all through 
galaxies at high redshifts. 

Previous applications cover the range from star clusters, 
normal galaxies E, ... Sd, dwarf galaxies with and with- 
out starbursts, tidal dwarf galaxies, interacting and merging 
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galaxies with their major starbursts, galaxy transformation 
processes in galaxy clusters, high redshift galaxies with and 
without starbursts and post-starbursts and damped Lyman- 
a absorbers. An early attempt at coupling GALEV evolution- 
ary synthesis models into a cosmo dynamical structure fo r- 
mation simulation was presented in lContardo et all (|l99&t l. 

GALEV models are now widely used throughout the com- 
munity. To facilitate access to the latest developments we 
here present a user- friendly and customized web-interface. 
It allows access to already available models for the evolu- 
tion of star clusters of various metallicities, and galaxies of 
all types both in terms of their time evolution for compar- 
ison with observations in the Local Universe and in terms 
of their redshift evolution. Furthermore it allows the user to 
run new models for specific applications. 

The philosophy for GALEV models is to keep them sim- 
ple with as small a number of free parameters as possible, 
and have them predict a large range of observational prop- 
erties, which - in comparison with observations - constrain 
the few free parameters. At the present stage, GALEV mod- 
els are 1-zone models without spatial resolution and without 
any dynamics included. Future prospects are the consistent 
inclusion of dust detailing absorption and reemission as a 
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function of gas content, metallicity and galaxy type, and the 
coupling with a dynamical model for stars and gas, includ- 
ing a star formation criterion and an appropriate feedback 
description, to cope with spatially resolved galaxy data. 



2 THE GALEV CODE: AN OVERVIEW 

2.1 Evolutionary synthesis for star clusters and 
galaxies 

Our GALEV evolutionary synthesis models have many prop- 
erties in common with the evolutionary synthesis codes 
from ot her groups, e.g. BC03 (iBruzual fc Charlothool ). PE- 
GASE (fFioc fc Rocca-Volmerangel Il997l ^ and Starburst99 
ijLeitherer et al.lll999L just to name a few, in that all these 
codes trace the evolution of the stellar population in terms of 
integrated spectra and/or colours for simple and composite 
stellar populations. 

In contrast to evolut ionary synth esis , stellar population 
synthesi s (e.g. [0'Connel|ll976l. Il980l) or differential syn the- 
sis (e.g. IPicklel l 1985al lbl; iPickles fc Visvanathanl I1985T ) at- 
tempts to find the best linear combination of stellar spectra 
from some library to fit an observed galaxy spectrum . This 
approach usually achieves very good fits but is limited to 
a status quo description and has difficulties to prove the 
uniqueness of its solutions. The existence of a stellar Initial 
Mass Function (IMF) and some continuous Star Formation 
History (SFH) can be imposed as boundary conditions via 
Lagrangian multipliers. The major advantages of the popu- 
lation synthesis approach are that it can give valuable first 
guesses for unknown SFHs and that it allows for unexpected 
solutions. 

All evolutionary synthesis models, on the other hand, 
have to assume a stellar IMF and a SFH, i.e. the time evolu- 
tion of the SFR for the galaxy. They use stellar evolutionary 
tracks or isochrones that have to be complete in terms of all 
relevant stellar evolutionary stages. 

Both methods need stellar spectral libraries that also 
have to be complete in terms of stellar effective temperatures 
T e ff, surface gravities log(g), and metallicities [Fe/H]. 

GALE V models a re ava il able for a range of stel lar IMFs, 
includ ing ISalpeter! (|l955T ), iKroupal <|200lh . and IChabrierl 
(|2003h . Other choices of the IMF can easily be customized. 

For the case of a simple stellar population (SSP), i.e. a 
star cluster, the SFH is a 5-function, meaning that all stars 
are formed in a single timestep. GALEV models for SSPs of 
different metallicities (iKurth et al.ll 19991; ISchulz et alj|2002l ; 



lAnders fc Fritzei 120031 ; iLillv fc Fritzei bOOd ) were shown to 
well reproduce observed colours and spectral indices of 
star clusters as function of age and metallicity. They also 
show a pronounced non-linearity at metallicities close to 
and above the solar value. Important results include the 
findings that any colour-to-age or index-to-age calibra- 
tion/transformations are only valid at one metallicity and 
that any colour-to-metallicity or index-to-metallicity cali- 
bration / transform ations are only valid at a given age (cf. 
ISchulz et al.ll20o3 ) , and that extrapolating observational re- 
lations beyond their calibration range lead to significantly 

mislea ding results. 

In lAnders fc Fritzei (|2003l ) we demonstrated the impor- 
tance of nebular emission lines and continuum for young 



stellar populations and showed that they can account for 
as much as 50 — 60% of the flux in broad-band filters, in 
particular at low metallicities. 

For the description of undisturbed galaxies, SFHs have 
been determined for normal average galaxies of types E, SO, 
Sa, Sb, Sc, Sd (cf. Sect. I5.2JI , that, in combination with a 
Salpeter IMF extending from a lower mass limit of 0.1 M@ 
(roughly the hydrogen burning limit) to an upper mass limit 
around 70 — 120 Mq , depending on the set of isochrones (cf. 
Sect. 13.1) selected, provide agreement with average observed 
galaxy properties in terms of colours, spectra, luminosities, 
abundances, and gas content. We stress that our galaxy 
types are meant to denote spectral types and we caution 
that the one-to-one correspondence between spectral and 
morphological types observed in the Local Universe might 
not hold to arbitrarily high redshifts. 

While in terms of spectral evolution of the integrated 
light of galaxies (or star clusters) our models are compara- 
ble to other evolutionary synthesis models, they go beyond 
those in that they also allow to describe and analyze re- 
solved stellar populations in terms of colour-magnitude di- 
agrams (CMDs) and in that they self-consistently describe 
the chemical evolution of the ISM in galaxies together with 
the spectral evolution of the stellar populations (for the lat- 
ter see Sect. 12, 3[) , allowing to realistically account for the 
coexistence of stellar subpopulations of different metallici- 
ties observed in local galaxies. 

Most of the aforementioned capabilities are not en- 
tirely new. Models of the photometric e volution of galaxies 
date back to iTinslevi (|l967l , Il968l . I1972T ). the fi rst spectro- 
scopi c models appeared roughly a decade later (|Bruzual A.I 
1 19831 ; iGuiderdoni fc Rocca-Volmerangel Il987fl. Models of 



the chemical evolu t ion of galaxies (e.g., iTruran fc Cameronl 
Il97ll ; lTinslevl[l972l ; iMatteucci fc Padovanilll993l 'l were capa- 
ble of taking the increasing enrichment of subsequent stel- 
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and line indices (e.g.. Weiss et al.lll995l;lBressan et al.lll996l). 


ISilva et alj (|l998h were the first to include a radiative trans- 



fer code into their model and could hence extend the wave- 
length coverage into the (far-)infrared. However, only very 
few models (e.g.. iPrantzos et al. l ll994l ; |Portinari et al.|[l998f ) 
exist that take the metallicity dependence of stellar yields 
into account and hence merit to be called chemically consis- 
tent. Those, unlike GALEV presented here, mainly focus on 
the metallicity distribution in the solar neighbourhood and 
only derive a detailed chemical evolution but no spectral or 
photometric evolution. 

2.2 Chemical evolution of galaxies 

Modeling the chemical evolution of galaxies starts from a 
gas cloud with given initial (e.g. primordial) abundances 
and given ma ss. A modified version of Tinsley's equations 
(|Tinslevlll968t ). including detailed stellar yields, is solved to 
study the chemical enrichment history of galaxies of different 
spectral types. This requires knowledge of stellar yields, i.e. 
production rates of different elements and isotopes, includ- 
ing contributions from SN la, as well as stellar lifetimes as 
function of stellar mass and metallicity, that can be taken 
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from nucleosynthesis and stellar evolution models, respec- 
tively. 

Closed-box models can be compared to models with 
specified in- and outflow rates and abundances. We follow 
the chemical evolution of a large number of chemical ele- 
ments H, He, ... Fe, fully accounting for the time delay be- 
tween SF and the return of material in stellar winds, PNe, 
and SNe. 



2.3 Chemically consistent GALEV models for 
galaxies 

Combining the chemical evolution of ISM abundances and 
the spectral evolution of the stellar population thus allows 
for what we call a chemically consistent treatment of both 
the chemical evolution of the ISM and the spectral evolution 
of the stellar population in galaxies: we use input physics 
(stellar evolutionary tracks, model atmospheres, stellar life- 
times and yields) for a large range of metallicities and consis- 
tently account for the increasing initial abundances of suc- 
cessive stellar generations. 

Broad stellar metallicity distri butions have been re- 
porte d for t he Milky Way disk iR ocha-Pint o fe Maciell 
119981 ) . bulge (ISadler et al.1 1 19961 ; iRami'rez et alJl200Ch . and 
halo (|Ak et al.l l2007tl . as well as for the nearby elliptica l 



galaxy NGC 5128 (|Harris et al.ll 19991 ; lHarris fe Harr5l200d) . 

Depending on the SFH of the respective galaxy type 
and eventually its infall rate, stars of different ages within 
a galaxy will have different metallicities and obey an age- 
metallicity relation determined by their galaxy's SFH. 

An important consequence of this coexistence of stars 
with different ages and metallicities is that stars of different 
metallicities and different ages dominate the light in different 
wavelength regions. It has severe implications for metallic- 
ity indicators defined in different wavelength regimes, which 
cannot be expected to trace one and the same stellar metal- 
licity. It also affects some widely used SF indicators and 
modifies, e.g., the calibrations for SFRs fr om Ha or [OH] 
fluxe s, as well as from FUV luminosities (cf. iBicker fe Fritzej 
120051 '). 

Our GALEV code can model the spectral and chemical 
evolution of galaxies with arbitrary IMF and SFH over cos- 
mological timescales, from the very onset of SF all through a 
Hubble time. In combination with a cosmological model we 
can follow the redshif t evolution of gala xies from the early 
universe until today (|Bicker et al.|[20o3 ). It also allows to 
directly study the impact of evolutionary corrections as well 
as of the chemically consistent treatment a s compared to us- 
ing so lar metallicity input physics only (cf. iKotulla fe Fritzej 
120091 ). 



2.4 Colour magnitude diagrams 

Despite simplifications the CMDs are valuable tools to 
study systematic effe cts o f SFH recovery from observations. 
iFritze fe Liilvl (|2007l ) and lLillv fe Fritzej (|200Sf ) investigated 
the accuracy of recovering SFHs from CMDs as a function 
of the ages of various subpopulations. In addition, they com- 
pared their results with the accuracies of recovering SFHs 
from integrated spectra, multi-band photometry, or Lick in- 
dices. These systematic studies are essential for more dis- 
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Figure 1. Stclar evolutionary tracks of a 2Mq star for different 
metallicities ranging between Z = 0.0004 = 1/50 Zq and Z = 
0.03 = 1.5 Z e . 



tant unresolved stellar populations. Moreover, since syn- 
thetic CMDs can be calculated in any desired filter combina- 
tions, they can be used to optimize observational strategies 
with respect to the optimal filter combination, e.g. to disen- 
tangle ages and metallicities of young, intermediate-age or 
old stellar populations. 



3 INPUT PHYSICS 

In the following sections we will review the input physics we 
use for our GALEV models. 



3.1 Stellar evolutionary tracks and/or isochrones 

Data for stellar evolution can be taken either from isochrones 
or stellar evolutionary tracks, both having their advantages 
and disadvantages. GALEV models currently use the most re- 
cent co nsistent set of theore tical isochrones from the Padova 
group (|Bertelli et al.l Il994j . ff ) for five different metallici- 
ties [Fe/H] = (-1.7; -0.7; -0.4; 0.0; +0.4) and include the 
TP-AGB phase, t he importance of which was shown in 
ISchulz et al.l (|2002l 'l. In order to be able to fully account for 
emission lines we also include the Zero Age Main Sequence 
(ZAMS) into our models. The isochrone for this ZAMS in- 
cludes stars up to 120 Mq and is created from the unevolved 
first data points of the stellar evolutionary tracks. F or more 
details we refer the reader to lBicker fe Fritz el (|2005l). 

In Fig. Q] we show stellar evolution tracks for a 2 Mq 
star for 6 different metallicities r anging from Z = 0.0 004 = 
1/50 Z Q to Z = 0.03 = 1.5 Z© (|Girardi et al.ll2000l ). The 
general shape remains unchanged with changing metallic- 
ity, but low-metallicity tracks are shifted towards increasing 
luminosities and higher effective Temperatures, i.e. towards 
the top left in the Hertzsprung-Russell diagram. This be- 
haviour is the same for stars of all masses, in the sense that 
with decreasing metallicity stars become more luminous and 
their spectra shift towards higher effective temperatures. It 
is therefore crucial to include those effects to obtain a con- 
sistent picture of galaxy evolution. 
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Figure 2. Stellar spectra from the Lcjeunc library for a star with 
identical T c fj = 3000K and logg = 4.0, but different mctallicities 
from [M/H] = -1.5 . . . + 0.3. 



3.2 Library of stellar spectra 

In principle, every library of stellar spectra - observed or 
theoretical - can be used, provided it is complete in terms 
of stellar T e g, log(g), and [Fe/H]. GALEV assigns spectra 
to stars in the isochrones according to each star's metal- 
licity, effective temperature T c g and surfce gravity log(g), 
normalising the spectra with the star's luminosity (for de- 
tails see Sect. 14. I'll . So far , GALEV uses the BaSeL lib rary of 
model atmospheres from iLeieune et al ] l| 19971 . 1 19981 ). origi- 
nally based on the iKurucj (|l992h library. The wavelength 
coverage spans the range from the XUV at A ~ 90 A to the 
FIR at A = 160 fim, with a spectral resolution of 20 A in 
the UV-optical and 50 — 100 A in the NIR range. We remind 
the reader that there are significant contributors other than 
starlight (e.g. PAHs, thermal emission from cold dust) at 
wavelengths beyond the K-band, that are not currently in- 
cluded in our models. This wavelength range should hence 
be used with caution. 

Stellar spectra are heavily influenced by metallicity due 
to the increased absorption line strength (e.g. for Fe-lines) 
and line-blanketing with increasing abundances. In Fig. [2] 
we show the example of a cool star (T c g = 3000K, logg = 
4.0) from the above mentioned BaSeL library. The effects of 
metallicity are stronger for cool stellar atmospheres where 
molecular absorption by VO, TiO, NH4, H2O etc. plays a 
larger role (jAllard et al 1 l2000l : iKucinskas et al.ll2006l ). 



3.3 Gaseous emission: lines and continuum 

In addition to the stellar absorption spectra we also compute 
line and continuum emission from gas ionized by hot massive 
stars. 

We do not take the ionizing photons from the stellar 
spectral library, but i nstea d use the tabulated values from 
ISchaerer fe de Koterl (|l997l ) for N Ly c as a function of stel- 
lar effective temperature, radius, and metallicity. They yield 
much better agreement with observations as well as with re- 
cent results from expanding non-LTE, line-blan keted models 
atmosphere calculations bv lSmith et all (|2002l) . 



We can then compute the total flux emitted by the gas 
per unit wavelength as 
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with the speed of light c, electron temperature 
T = 10 000 K, and total recom bination c oefficient 
a (2) (10 000K) = 2.575 xl0~ 13 cm 3 s _1 (jAllerlll984l ). The gas 
continuu m coefficients are th e n com puted for T = 10 000 K 
following lErcolano fc Storevl (|2006l ) that contains an algo- 
rithm to compute the bound-free radiation for hydrogen 
and helium (7m, THei, and 7Heii). The HI two-photon emis- 
sion c oefficient 7h2 P is taken from iNussbaumer fc Schmutzl 
(Il984j) and for the free- f r ee em ission 7b we use the formula 
from Brown fc Mathews] (Il970j) and compute the Gaunt fac- 
tors using the algorithm from Irlummerl (|l98ST ). All 7- factors 
are summed to form the final: 
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7b + 7H2p + 7HI 



n Hc + 



-7HcI 



n He+ + 



-7HcII- 



(2) 



n H + n H+ 

For the densities of the helium He I and He II ions rela- 
tive to HI we use values typical for HII regions in galaxies 
(Ercolano, priv. comm.), 0.0897 and 1.667 x 10" 4 for He I 
and Hell, respectively. The final isochrone spectra are not 
very sensitive to slight changes in these values. 

The line strengths of the hydrogen lines are computed 
using atomic physics and the emission rate of ionizing pho- 
tons NLyC from O- and early B-stars. From the number of 
ionizing photons, NLyC, we compute the flux in the H/3 line 
using 



f(H/3) = 4.757 x 10" 13 erg x N LyC . 



(3) 



Assuming Case B rec ombination for a pure hydro gen cloud 
with T c = 10 000 K (|Osterbrock fc Ferlandl 120061 ). we can 
then derive the line strengths of all other hydrogen lines. 
Line strengths for heavier elements are computed using 
metallicity-dependent line-ratios relative to H/3. For metal- 
licities [Fe/H] > —0.4 those are taken from the IStasinskal 
(| 19841 ) photoionization models, adopting typical Galac- 
tic values of T c = 8100 K and n e = 1 cm -3 for electron 
temperature and density. For lower m e tallici ti es, w e use 
observed line ratios f rom llzotov et al.l (| 19941 . Il997l ) and 
llzotov fc Thuanl (|l998l ) . These are supposed to include sys- 
tematic changes in T c , n e and the ionizi ng radiation field 
in lo wer metallicities environments (cf. lAnders fc Fritze] 
I2003D . Note that we do not include detailed radiation trans- 
fer calculations (as . e.g.. in [Garcia Vargas fc Diaz! 1 19941 : 
iFerland et ail 1 19981 : lErcolano et al.l |2003j ). hence we cau- 
tion the user to use line-ratios for detailed line diagnostics. 
We correct the gaseous emission for small amounts of dust 
within the HII regions by reducing the ionizing flux by 30% if 
the gaseous metallicity is near solar, i.e. for [Fe/H] ^ —0.4. 
For lower metallicities we use the full ioni zing flux sinc e 
those environments are essentially dust-free l|Mezgerlll978l ). 
We do not account for yet unknown amounts of dust deple- 
tion of heavy elements within the HII regions. However, in 
the case of low-metallicity galaxies these effects are already 
included in the observed line-ratios. 

Note that we do not account for internal self-absorption 
of Lyman continuum photons within the HII regions or 
the surrounding galaxy, since the fraction of flux leak- 
ing out of th ese regions is still a matter of ongoing de- 
bate (see, e.g. jFerguson et al.ll 19961 : ICastellanos et al.ll2002l : 
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Figure 3. Emission line spectrum of the 4 Myr isochrone for 5 dif- 
ferent metallicities ranging between Z = 0.0004 ([Fe/H] = —1.7) 
and Z = 0.05 ([Fe/H] = +0.4). 



Fernandez-Soto et all 120031 : llnoue et all 120051 : ISiana et all 
20071 ; IWise fc Cenll2008l ) 



Depending on application, our models can include both 
continuum and line emission, continuum emission only 
or no gas emis s ion a t all. In iKriiger et alj l|l995l ) and 
lAnders fc Fritzd (120031 ) we showed that in the case of Blue 
Compact Dwarf galaxies and in young and metal-poor SSPs, 
gaseous emission can contribute as much as 60% to the flux 
in broad-band filters. Emission lines are the dominant con- 
tributors in the optical, whereas continuum emission domi- 
nates in the NIR. 

Combining the effects of higher luminosities and higher 
effective temperatures with the effects of longer lifetimes of 
high-mass stars at low metallicities has a profound impact 
on the spectrum of galaxies. All three factors lead to bluer 
colours, higher overall luminosities and as a further aspect 
significantly stronger gaseous emission, i.e. gaseous contin- 
uum and line emission. 

In Fig [3] we show the gaseous emission spectra for 
isochrone spectra at an age of 4 Myr and for 5 metallicities 
from 1/50 Zq to 2.5 Zq. The first notable aspect is the over- 
all emission line strength that is higher by a factor of about 
10 comparing the two extreme cases. But also the individ- 
ual line ratios change significantly, i.e. the ratio between Ha 
(6563 A) and [N II] (6583 A) (Ha/ [Nil] « 48 at low metal- 
licity and ~ 7 for solar metallicities). This extreme ratio 
is directly affected by the lower nitrogen abundance at low 
metallicity. Other line ratios, e.g. [O II] (3727 A) to [O III] 
(5007 A) also change due to the more intense radiation field 
coming with the low metallicity environment. 



3.4 Lick stellar absorption indices 

Since the resolution of the iLeieune et al.l (| 19971 . Il998l ) li- 
brary is not sufficient to calculate Lick absorption indices di- 
rectly from the spectra, GALEV mode l s use t h e empirical cal- 
ibratio ns presented by Gorgas et al. | l993h , IWorthev et al.l 
(|1994| ) and lWorthev fc Ottavianil l| 19971 ) .Rom those, GALEV 



calculates the fluxes in the Lick indices for every star, sums 
them up for the entire stellar population at every timestep 
to yield integrated index fluxes and, in combination with 
continuum fluxes from the integrated absorpti on line spec- 
tra, t he respective index equivalent widths (|Kurth et al.l 
Il999l : iLillv fc Fritzd l200rj| ). GALEV at the present stage in- 
cludes the following Lick indices: H<5a, R~fA, H<5f, Hy.F, CNi, 
CN 2 , Ca4227, G4300, Fe4383, Ca4455, Fe4531, Fe4668, H/3, 
Fe5015, Mgi, Mg 2 , Mgb, Fe5270, Fe5335, Fe4506, Fe 5709, 
Fe5782, NaD, TiOi, and Ti0 2 (see lTrager et alJll99Sl . and 
references therein for all index definitions). 



3.5 Stellar yields 

To model the chemical enrichment histories of galaxies, 
GALEV uses stellar yields for a large number of individ- 
ual elements (H, He, Li, Be, B, C, N, O, F, Ne, Na, Mg, 
Al, Si, P, S, CI, Ar, K Ca, Sc, Ti,V, Cr, Mn, F e, Co, Ni, 
Cu, Zn, Ga, Ge) from IWooslev fc Weaver! (Il995l) for mas- 
sive stars and from lvan den Hoek fc Groenewegen (|l997l ) for 
low-mass stars of various metallicities. Stellar lifetimes are 
taken from the isochrones. There is a min or inconsistency in 
doing so, since IWooslev fc Weaver! (| 19951 ) used models with- 
out mass-loss, while current isochrones in general account 
for mass-loss. However, since yields are only available for a 
very coarse metallicity grid this does not significantly affect 
the resulting chemical evolution. GALEV also includes type 
la SN yields for th e carbon deflag r ation whit e dwarf binary 
scenario (W7) from lNomoto et all (|l997l ). See lLindner et al.l 
( 1999) for a detailed description of the chemically consistent 
chemical evolution aspects of GALEV. 



3.6 Ejection rates and remnant masses 

One of the central input parameters for GALEV is the time- 
dependent ejection rates necessary to compute the chemical 
evolution of galaxies. Those rates are derived from the initial 
stellar masses M* and the remaining remnant masses hir. 
For stars with masses M* ^ 30 M the remnant is assumed 
to be a black hole of meH = 8.0 M©, with the remaining 
mass being returned to the ISM; stars with initial masses of 
30 Mq ^ M* ^6.0 Mq result in a neu tron star of mass uins 
given by iNomoto fc Hashimoto! (|l988l ): 

m NS [M Q ] = 1.02 + 3.6363 x 10~ 2 (M*/M Q - 8.0)M Q ; 

For the mass range of 6.0 Mq m* ^ 0.5 Mq, for which the 
stellar re mnant is a white dwarf, we use a fit to the data 
points of lWeidemannl l|2000l ): 

mwD [Mq] = 0.444 + O.O838(M*/M ); 

This, combined with the extrapolation of the NS relation 
down to 6Mq, provides a better matched connection be- 
tween the two mass ranges, whil e being compa t ible w ith the 
slightly steeper slope derived by iKalirai et al.l (|2008l ). How- 
ever, as the remnant is only used to derive the mass returned 
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Figure 4. Main-sequence lifetimes as function of initial stellar 
mass for three different mctallicities. 



to the ISM during its lifetime and death, the exact transition 
point from neutron star to white dwarf is of minor impor- 
tance. 

Stars with masses m* < 0.5Mq have lifetimes in excess 
of a Hubble time and negligible winds, hence do not return 
any material to the ISM. 

In Fig. U we show the main-sequence lifetimes of stars 
as function of initial mass for three different metallicities 
from Z = 0.002 = 1/10 Z s to Z = 0.02 = Z , based on data 
fromlMarigo et all (|200St ) for M t < 7M and lBertelli et all 
(| 19941 ) for M+ > 7Mm. Low-mass stars live longer at high 
metallicities, while high-mass stars have longer lifetimes for 
lower metallicities; stars with masses ~ 3 — 4Mq have 
roughly the same lifetimes independent of metallicity (see 
Fig. [31 . The 20Mq star for example has a 10% increased 
lifetime at Z — 1/10Zq as compared to Z = Z©. 



3.7 Filter functions and magnitude systems 

GALEV includes a large number of filter functions to be con- 
volved with the model spectra in order to avoid uncertain- 
ties associated with transformations between filter systems. 
It also provides the option to choose the desired magni- 
tude system Vegamag, ABmag, and STmag to be directly 
comparable to observations and avoid transformations be- 
tween different magnitude systems. Magnitudes in the Vega- 
mag systems are defined to have magnitud e zero for an 
AOV s t ar; we use the Vega-spectrum from the lLeieune et ail 



|l997 , 19981) library combined with the flux calibration 



from 



Bohlin fc Gillilandl (l2004h as our standard star. AB 



magnitudes |Okel Il974 iBohlin fc Gillilandl 120041 ) are de 



rived from the monochromatic flux f v such that juab = 
— 2.51og(f„) +48.6 if f v is measured in erg s _1 cm -2 Hz -1 ; 
a colour of in the AB magnitude means that the object 
emits constant flux per unit frequency interval; analogous 
to the above, colours of in the ST magnitude systems 
mean constant emitted flux per unit wavelength interval. 
The zeropoint has been chosen such that a source with 
fA = 3.63 x 10~ 9 erg s" 1 cm" 1 A -1 has hist = Omag in all 
filters. 

Our current database contains filters from all major 



telescopes and instruments, including all HST instruments 
(incl. WFC3), many ESO instruments, and all common fil- 
ter sets like Johnson/Cousins, Stromgren, Washington, and 
SDSS. However, note that due to the wide wavelength sam- 
pling of 20 A of the current spectral library, narrow-band 
filters cannot be adequately supported at this stage. 



3.8 Dust extinction 

To account for extinction and reddening due to interstel- 
lar dust, GALEV also implement s the most commonly used 
empirical extinction laws from ICalzetti et al.l (|2000f l and 
ICardelli et all (|l989h . The earlier was derived from actively 
star-forming galaxies and describes a relatively gray extinc- 
tion without the 2175 A bump, characteristic for the latter 
extinction law. For the Cardelli law we assume a mean ex- 
tinction parameter of Rv = A(V)/E(B — V) = 3.1, char- 
acteristic for relatively quiescent galaxies. Although models 
exists that offer a more detailed treatment of dust extinc- 
tion and in some cases dust emission ( s ee, e. g., Silva et al 
ll99Sl ; |Popescu et alJIiuOOl; ICunowll200ll. bOO^lPopita et al 



20051 : (Mollenhoff et al.ll2006l ; IPiovan et al.ll2006al lbl) . the 



jor drawback of these models is that they in general require 
assumptions on the geometry and spatial distribution of 
dust, gas and stars. However, the aim of GALEV is to de- 
scribe the spectrum of the average represantation of each 
galaxy type, making if difficult to compare these two ap- 
proaches. 



3.9 Cosmological model 

GALEV can also be coupled to a cosmological model to de- 
scribe the evolution of galaxies as function of redshift. For 
this purpose we have to convert ages into redshifts and vice 
versa. In its current version we implemented a flat cosmology 
(SIk = 0). The choice of the local Hubble constant Ho, the 
density parameters Qm and J1a (with the additional con- 
straint Am + = 1) and the galaxy formation redshift 
Zform then completely determine the galaxy age as function 
of redshift. Because of the short time interval between e.g. 
z = 12 and z = 5, the exact value of zf orm has very little im- 
pact, hence we choose an intermediate value of Zf orm — 8 as 
our default. 



3.10 Intergalactic attenuation 

To correctly describe spectra of higher redshift galaxies one 
has to include the absorption shortwards of the Lyman-a 
line due to intervening neutral hydrogen clouds. For that 
reason GALEV implements the description for the average 
attenu ation effect as function of redshift following iMadaul 
|l995l ) and covering the range ^ z 7. 



4 PROGRAM STRUCTURE 

The actual modeling process with GALEV can be divided into 
three steps: 

(i) In the first step GALEV convolves the isochrones of 
each age- metallicity combination with the specified IMF, 
normalized to 1 M© . It then assigns a spectrum from the 
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stellar library to each star on each isochrone and computes 
the integrated isochrone spectrum. Then the gaseous con- 
tinuum emission and emission lines are added to the young 
isochrone spectra. Using the yield tables and stellar remnant 
masses described in Sect. 13.51 with stellar lifetimes from the 
isochrones, GALEV also derives the gas and metal ejection 
rates needed for the chemical evolution. 

(ii) In the second step GALEV computes the chemical and 
spectral evolution of the desired stellar population. For each 
timestep the current isochrone spectra are weighted with the 
SFH. The contributions from older SF episodes are obtained 
by integration over all past timesteps. The required interpo- 
lations in age and metallicity are described below. GALEV 
thus calculates the time evolution of the integrated spec- 
trum and its resulting line strength (i.e. Lick-) indices of a 
simple (SSP, star cluster) or composite (galaxy) stellar pop- 
ulation, including the gaseous emission where appropriate. 
From the ejection rates and the available gas mass GALEV 
calculates the new gaseous metallicity that will be used for 
stars born during the next timestep. 

(iii) In the last step GALEV converts the integrated spec- 
tra into magnitudes in a large number of filters. Given a list 
of requested filters it convolves the spectra with the filter 
functions and applies zero-points to yield absolute magni- 
tudes for each timestep. It can alternatively be combined 
with a cosmological model to yield apparent and absolute 
magnitudes as function of redshift. If requested it also ac- 
counts for dust extinction using observed extinction laws 
and, in the context of cosmological evolution, also for the 
attenuation by intergalactic neutral hydrogen. 

This process is also explained in a more vivid step-by-step 
explanation in appendix 1X1 



4.1 Computation of isochrone spectra 

A crucial step is the assignment of stellar spectra from the 
library to the points describing the isochrones. Since the 
available stellar parameters T e ff and log(g) in the library 
often do not match those required by the isochrones, GALEV 
has to interpolate between them in both T e ff and log(g). 
For a given combination of T c fr and log(g) this is done as 
follows: 

(i) Find up to four spectra bracketing the required values 
in both T e ff and log(g). 

(ii) Interpolate the spectra to the required value of T e ff 
in each pair of lower (upper) values of log(g), yielding two 
new interpolated spectra with the correct T e ff and different 
log(g). An important factor during this interpolation is the 
weighting of the spectra with each star's luminosity given 
by the isochrone. We choose to use the integrated luminos- 
ity in the Johnson-V or Bessell-H-band, depending on the 
temperature of the star given by T c ff. The original approach 
to normalize all stars in the V-band turned out to be insuf- 
ficient for cool giants with only little flux in the optical. A 
"cool giant" in this context is defined by T c a ^ 3500 K and 
log(g) < 3.5. 

(iii) The spectrum for the required value of log(g) is then 
obtained by interpolation between the two spectra with the 
right T c ff, again weighting with each star's respective lumi- 
nosity. 



For the very hot stars (T e ft > 50 000 K), the BaSeL stellar li- 
brary does not provide spectra. In those cases we extend the 
spectral library by approximating the missing stellar spectra 
with black-body spectra of the requested temperatures. The 
validity of this approximation for wavelengths longwards of 
A ~ 230 A is supported by only minor differences between 
pure black-body a nd true spectra from both o bservations 
jGauba et al.ll200lh and modelling (|Rauchll2003l ) of very hot 
central stars of planetary nebulae. 

In a final step, all isochrone spectra are normalized to a 
distance of 10 pc, and are given in units of erg s _1 cm~ 2 A -1 . 

4.2 Interpolation of isochrone-spectra and 
integration of galaxy spectra 

One important aspect of GALEV is how to interpolate be- 
tween the isochrone spectra of different ages and metallici- 
ties, in other words how to map the coarse grid of isochrones 
available onto the finer grid needed for galaxies. GALEV here 
interpolates logarithmically between ages and linearly be- 
tween metallicities expressed as \og(Z) ~ [Fe/H]. During 
their very early stages the gaseous metallicities of galax- 
ies are lower than the metallicity of the lowest metallicity 
isochrone with [Fe/H] = —1.7. For stars born at these stages, 
we chose to use the lowest metallicity isochrone to avoid 
the uncertainties in an extrapolation to lower metallicities. 
The same is done for stars born from gas with a metallic- 
ity higher than the highest metallicity isochrone. For these 
stars we use the highest metallicity isochrone available with 
[Fe/H] = +0.4. 



5 CALIBRATION OF THE GALEV MODELS 
AND COMPARISON TO OBSERVATIONS 

We stress that with the input physics as outlined above, a 
stellar IMF with lower and upper mass limits chosen and 
total mass of the galaxy or star cluster specified, GALEV 
models provide the time evolution of spectra, luminosities, 
and colours in absolute terms. The same holds true for the 
gas content and the chemical enrichment of galaxies. No a 
posteriori gauging is applied. 

The only exception is the cosmological context for 
galaxies: before luminosities and magnitudes of redshifted 
galaxies are calculated, the B-band model luminosities af- 
ter a galaxy age corresponding to redshift z = 0.0044 (i.e. 
the redshift of the Virgo cluster) in the chosen cosmological 
model are scaled to the average observed B-band luminosity 
of the respective galaxy type in the Virgo cluster (cf. Sect. 
E3J. 

In the following we present as an overview the compari- 
son of GALEV models for star clusters and for normal galaxy 
types E, Sa . . . Sd with observations and refer to previous 
papers for more details. 

5.1 Star clusters 

The colours from U through K predicted by GALEV mod- 
els for SSPs or star clusters at an age of ~ 12 — 13 Gyr 
are in very good agreement with the respective observed 
colo urs for a large set of M31 and Milky Way glob ular clus- 
ters (jBarmbv et alj|200d : iBarmbv fc Huchrall200Ch at their 
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respective metallicities (cf. ISchulz et al.ll2002l ). Small devia- 
tions in the U- and B-bands can be explained by the ex- 
istence of Blue Straggler stars (likely products of stellar 
mergers) in the dense cores of GCs, which can make sig- 
nificant contributions at those short waveleng ths (|Xin et al.l 
120071 : IXin fc Penal 120051 : ICenarro et all 120081 ). but are not 
included in our models yet, as standard isochrones exist 
for non-interacting single stars only. Further contaminants 
in the blue-to-UV reg i on are Blue Horizo ntal Branch stars 
jSchiavon et all |2004| ; iDotter et all 120071 ) and to a lesser 
degree even Cataclysmic Va riables and low-mass X-ray bi- 
naries in the far-UV range ijRich et all 1 19931 iDieball et all 
l2007h . 

As shown in iKurth et all (|l999h and iLillv fc Fritzel 
(|2006l ). Lick indices for SSPs as calculated by galev agree 
well with Galactic an d M31 globular cluster data as com- 
piled by (|Harrislll996f ) in his online catalog. 

GALEV models for star clusters also show good agree- 
ment at an evolutionary age of ~ 12 — 13 Gyr with the em- 
pirical calibrations for colours (B — V) and (V — I) versus 
metallicity [Fe/H] f or old Galactic and M31 globular clus- 
ters, a s e.g. given bv lCouture et all (|l990h and lBarmbv et all 
(|2000T ) over the metallicity range -2.3 ^ [Fe/H] < -0.5 of 
these clusters. They also show, however, significant devia- 
tions from a linear extrapolation of these empirical relations 
towards higher metallicities, and they also show that the em- 
pirical relations are only valid for the old globular clusters 
for which they have been derived. Models can, of course, be 
used to study the behaviour of colour - metalli city relations 
for a ny colour and as a function of time (cf. ISchulz et all 
12001 for details). 

In Figure [5] we plot the evolution of three SSP models 
with different metallicities ranging from the lowest available 
value of [Fe/H] = -1.7 to the solar value [Fe /H] = 0.0. For 
compa rison we also show V-K colours from iPersson et all 
(l983J) for star clusters in the Large Magellanic Cloud, 
colour- coded according to t heir m etallicity as derived from 
CMDs iMackev fc Gilmorel l|2003h . Our models are able to 
not only reproduce the full range of observed colours, but 
also match the colour-evolution of each metallicity subpop- 
ulation if accounted for typical uncertainties of 0.4 dex in 
the age-determination. For the old and metal-poor globular 
clusters our predicted V-K colours are in good agreement 
with observations. 

As shown in ILillv fc Fritzel (|2006h . galev models for 
SSPs also agree, after 12 — 13 Gyr of evolution with empirical 
calibrations of Lick indices vs. [Fe/H]. However, they also 
show that all metal-sensitive indices are also age-dependent 
and that the famous age-sensitive Hp and H-y indices are 
also metal-dependent t o some extent (see also lThomas et all 
120031 : iKorn et al.ll2005r i . Hence, the empirical calibrations of 
Mg2, Mgb, etc. vs. [Fe/H] are valid only for old (i.e. > 10 
Gyr) globular clusters. 

Empirical calibrations for colours or Lick indices vs. 
metallicity or age should therefore not be used for star clus- 
ters or globular clusters for which it is not a priori clear 
that their properties fall within the range of the calibrating 
Galactic clusters. Instead a full set of colours and/or indices 
in comparison to an extended grid of models covering the 
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Figure 5. Evolution of the V-K colour of our SSP models with 
metallicities of [Fe/H] = -1.7 (red, solid line), [Fe/H] = -0.7 
(green, dashed), and [F e/H] = 0.0 (blue, dash-do tted). We also 
plot the data ( ages f rom lMackev &; G ilmorc (20ojl), colours from 
IPersson et all l ll983h ) for 35 clusters in the Large Magellanic 
Cloud for comparison. 
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Figure 6. Star formation histories for the five different galaxy 
types E and Sa through Sd. Absolute numbers are normalized to 
a total galaxy mass of 10 10 Mq. 



full parameter space, allows for independent and accurate 
determinations of both ages and metallicities (cf. Sect. I8.1]) . 



5.2 Galaxies 

By de fault, all our g alev models for galaxies use a Salpeter 
IMF (|Salpeterlll955l ) with lower and upper mass limits of 
0.1 M Q and 100 M , respectively. 

The star formation histories, i.e. the time evolution of 
the star formation rates, of the different spectral galaxy 
types E and Sa to Sd are the basic parameters of GALEV 
and, in fact, of any kind of evolutionary synthesis models. 
In Fig. [6] we show SFHs for galaxies of different types from 
E through Sd, assuming the same total (i.e. stars and gas) 
mass Mtot of 10 10 M for all of them. These SFHs are in 
good agre ement wi t h chem ical and spectrophotometric find- 
ings from ISandagel (|l98rj ) . 

For the elliptical model we use an exponentially declin- 
ing SFR 
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*(t) = — x exp(-t/r) 



(4) 



with an e-folding time r = 1 Gyr and a = 8.55 x 10 yr. 
For spectral types SO and Sa through Sc the SFR is tied 
to the evolving gas content as given by 



*(t) = X 



M gas (t) 
10 9 : 



(5) 



with efficiency parameters p decreasing from early to 
later spiral types, P = 1.0 yr -1 for SO, P = 0.33 yr -1 for Sa, 
P = 0.19 yr -1 for Sb, and P = 0.11 yr^ 1 for Sc galaxies. 

Sd model galaxies are described by a constant SFR: 

M to 



*(t) = ip x 



10 10 



const, 



(6) 



where * = 0.4 M Q yr -1 . 

For closed-box models, these parameters a, P, ipo, and r 
are the only free parameters. We do not require additional 
parameters such as infall of gas or outflow in galactic winds 
to reproduce observations. Note that all these SFH param- 
eters are independent of galaxy mass. We therefore do not 
reproduce mass-metallicity or colour-magnitude relations for 
galaxies of identical spectral types. 

These SFHs are very si milar in all evo l ution- 
ary synthesis models (cf . e.g. iBruzual fc Chariot] 120031; 
iFioc fe Rocca-Volmerangel Il997l ). In detail, we adjust the 
SFH for the chemically consistent GALEV models as to 
match, after a Hubble time of evolution (i.e. « 13 Gyr 
for our assumed concordance cosmology with Ho = 
70kms~ 1 Mpc -1 , Q M = 0.30 and fl A = 0.70), the observed 

• average integrated colours from UV through NIR, 

• average gas fractions, 

• average metallicities, 

• average present-day SFRs, 

• average mass-to-light ratios, and 

• template UV - optical spectra 

of the respective spectral types as detailed below. All these 
observational constraints together very neatly define the av- 
erage SFHs of undisturbed galaxies E, Sa, Sb, Sc, Sd and 
tightly constrain the few parameters describing them for a 
given IMF. 

We stress that GALEV and all other evolutionary syn- 
thesis models with their respective SFHs are meant to de- 
scribe spectral types of galaxies. And we caution that while 
in the Local Universe and for undisturbed galaxies a clear 
one-to-one correspondence is observed between spectral and 
morphological types, it is an open question how far back in 
time this correspondence might hold. 

Gas fraction: The above described parameters have 
been tuned to reprod uce the typical gas fract ions observed 
in local galaxies (e.g. iRead fc Trenthamll2005l ). We use gas 
contents, defined as fractions of gas relative to the total 
(gas+stars) mass, of M gas /M t ot = 0.0 for E-models, 0.05 for 
Sa, 0.15 for Sb, 0.30 for Sc, and 0.55 for Sd models, respec- 
tively. 

Colours: With these parameters being fixed, we com- 
pare a wide range of model-predicted colours from near-UV 
to near-IR to values from the literature. Our (B-V) colours, 
e.g. of 0.86 (E), 0.78 (Sa), 0.64 (Sb), 0.56 (Sc), and 0.43 
(Sd) compare very well with the ranges found in the RC3 
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Figure 7. Evolution of galaxies of different types E through 
Sd in the (B-R)-(V-I) colour-colour-plane. Colours after 13 Gyr 
are plotted with large symbols. We also show colours for typical 

jalaxies (la rge circles around galaxy type) taken from lButa et all 

19951) and lButa fc William j Vl995f) . 



catalog of lde Vaucouleurs et al.l |l995T ). listing mean colours 



of (B - V) = 0.89^ (E), 0.74+°^ (Sa) 
0Mt°of (Sc), and 0.44±°;?* (Sd). 



0.66 



+0.22 



(Sb), 



Spectra: In Fig. 15.21 we compare our model spectra to 
local template s for the g a laxy t ypes E, Sa, Sb, Sc, Sd from 
the catalog of iKennicuttl (1 19921 ). As can be seen from those 
plots, the observed spectra for all spectral galaxy types from 
old passive ellipticals through the actively star forming Sd 
galaxies are well reproduced. Although our spectra (here us- 
ing Lejeune's library) have a lower spectral resolution than 
the observed templates, they nonetheless reproduce all fea- 
tures, like absorption and emission lines. In the lower part 
of each plot we also show relative differences between model 
and template spectra. Deviations are generally smaller than 
differences between different galaxies of the same type, con- 
firming the good agreement. Small differences for the emis- 
sion lines originate in different spectral resolutions of tem- 
plates and model spectra. Note that due to our chemically 
consistent treatment GALEV models can reproduce the Sd 
template spectrum with our Sd model at an age of 13 Gyr. 
This is a notable difference to all other evolutionary synthe- 
sis models that can only reproduce Sd te mplate spectra with 
youn ger models of ages 4 - 6 Gyr (cf. IBruzual fc Charlotl 
1 19931 ). 

Metallicities: As mentioned above GALEV models cal- 
culate the time evolution not only of individual element 
abundances but also of the global ISM metallicity Z. Af- 
ter ~ 13 Gyr of evolution (for d etails on the evolu tion 
as function of age or redshift see iBicker et al. | (|2004h or 
iKotulla fc Fritzd (|2009h ). our models reach ISM abundances 
of Z E = Z , Zsa = 1.5 Z , Zsb = 0.8 Z , Z Sc = 0.5 Z Q , 
and Zsd = 0.25 Z . Those metallicities are in good agree- 
ment with typical ISM abundances, as measured at 1R c j, 
the effective or c haracteristic radius, of spiral galaxies 
of various types (lOev fc Kennicuttl 1 19931 ; iFerguson et al.l 
1 19981 ; Ivan Zee et al.l 19981 ). Note that our models aim to 
represent the typical L* or M* galaxy of each spectral 
type. Galaxies of different masses within the same spec- 
tral type are known to have slightly different metallici- 
ties, as described by the mass-metallicity or luminosity- 
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Figure 8. Comparison of mod el spectra of diffe rent types (top row: E and SO, middle row Sa and Sb, bottom row Sc and Sd) with 
observed template spectra from iKenn icutt (1993). The black solid lines show the spectrum of GALEV models at an age of 12.8 Gyr, the 
coloured lines show spectra of several template galaxies of each respective type. In the lower part of each plot we also show relative 
differences between the GALEV spectra and each template spectrum. 
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meta l licity relations |s"killman et alj 1 19891 : iTremonti et al.l 
12004 : iKewlev fc Ellison! |200Sj ). As stated before, these re- 
lations are currently not accounted for by GALEV standard 
models. However, the user is free to choose parametriza- 
tions of the SFHs a s a function of galaxy mass (as, e.g., in 
iBressan et al.|[T994l ) that reproduce these relations. 

Star formation rates: GALEV models predict the fol- 
lowing present-day SFRs ty p in [Mq yr _1 ] for the different 
spectral galaxy types: ty p ~ for E, = 0.65 for Sa, 
* p = 0.4 for Sb, * p = 0.47 for Sc, and * p = 0.13 for Sd 
galaxies. These numbers are in agreement with the average 
star formation rates of 0.22, 0.32 ± 0.2, 0.33 ± 0.14, and 
0.09 ± 0.03M©yr _1 for Sa to Sd galaxies from iKennicuttl 
(| 19831 ). We caution that the smooth SFHs in GALEV and 
other evolutionary synthesis models certainly are simplifica- 
tions; real galaxies will in general experience fluctuations in 
their SFRs around these mean values. For the Milky Way, 
iRocha-Pinto et all (|200Cf) showed on the basis of individual 
stellar ages and abundances that the mean global SFH in- 
deed fell nicely between those of our Sb and Sc models, with 
fluctuations of a factor of 2 on timescales of 100 Myr. Such 
fluctuations around the mean, however, do not significantly 
affect the long-term evolution. 

Mass-to-light ratios: Using the stellar mass assem- 
bled after ~ 13 Gyr and the absolute V-band luminos- 
ity, we compute the V-band mass-to-light ratios. The re- 
sults are: M/L V (E, Sa, Sb, Sc, Sd) = 11.8,8.2,5.9,4.5,3.0, 
respectively. These include a factor that we call fraction 
of visible mass (FVM) and that accounts for the frac- 
tion of mass locked up in sub-stellar objects like brown 
dwarfs, planets, etc. This factor is set to FVM=0.5 in GALEV 
models, implying that only half the star formation rate 
forms luminous stars, while th e remaining half is locked 
up in non- luminous objects (cf. iBahcall et alj [l992h . Only 
with this F VM=0.5 can agreement with observed average 
M/L ratios dFaber fc Gallagher! 1 19791 : iBell fc de Jonell200ll ; 
iRead fc Trenthaml 120051 ) be achieved together with agree- 
ment in colours, spectra, and chemical abundances. Note 
that FVM does not include remnants of stellar evolution 
like white dwarfs, neutron stars or black holes; those are 
considered separately as part of stellar evolution. 

Introducing the FVM parameter does not only affect 
the mass-to- light ratios but also has an impact on the chem- 
ical evolution with in turn impacts on most other derived 
parameters. A Salpeter IMF with mass-limits of 0.1 and 
100 Mq would lock up 60% of the mass in long-lived stars 
with M ^ 1 Mq , with the remaining 40% progressing the 
chemical enrichment of the galaxy. However, by applying a 
FVM=0.5 we half this fraction of stars with M > 1 Mq so 
that now only 20% of stars return their ejecta at the end of 
their lifes, with the majority of mass being locked up. This 
concept to slow down the chemical enrichment is a common 
feature of many codes, but comes with diffe r ent n ames, e.g. 
as free parameter (1 — £) in lPortinari et ail (|l998l ). 

Luminosity normalization With the SFHs as given 
above and a Salpeter IMF, GALEV models calculate the time 
evolution of spectra and luminosities in absolute terms as a 
function of the total mass, initially all in gas. Colours, of 
course, are independent of galaxy mass as GALEV models 
assume the SFH of each galaxy type to be independent of 
galaxy mass (see Fig. [6j . Chemical ISM abundances are also 



determined by this SFH and by the absolute level of SFRs 
in relation to the mass in the initial gas reservoir. Spec- 
tral fluxes and luminosities on the other hand are mass- 
dependent. If models of undisturbed galaxies are to be com- 
pared to observations of individual galaxies, it is important 
to calibrate their absolute luminosities to the observed ones, 
thereby also calibrating their (gaseous + stellar) masses and 
their SFRs. 

Before we can compare models to distant galaxies, we 
need to calibrate their apparent magnitudes in the Johnson 
B-band against average observed magnitudes of the respec- 
tive spectral galaxy types in the Local Universe. The average 
observed apparent luminosities of the different ga l axy typ es 
in the Virgo cluster are given by ISandage et alj (|l985al lbl) 
to be m B (E) = 10.20, m B (S0) = 12.60, m B (Sa) = 11.95, 
m B (Sb) = 12.90, m B (Sc) = 12.83, and m B (Sd) = 13.95. The 
correction factors to match the computed and observed lu- 
minosities are then also applied to all mass-dependent model 
parameters, i.e. to stellar and gaseous masses, SFRs, etc. 



5.3 Disturbed galaxies: Starbursts &: SF 
truncation 

In addition to the undisturbed galaxies above with their 
smooth and monotonically decreasing SFRs, GALEV can also 
describe the evolution of galaxies encountering fast changes 
in their SFR, e.g. a starburst with significantly enhanced 
star formation rate, or a truncation or strangulation of the 
star formation rate on some shorter or longer timescale, re- 
spectively. This allows us to study not only the spectral 
but also the chemical evolution of galaxies that experience 
starbursts or a quenching of their normal SFR, as they can 
occur e.g. during galaxy-galaxy interactions or mergers or in 
the course of the various transformation processes discussed 
for field galaxies that fall into galaxy clusters. We recall that 
GALEV models only deal with effects on the SFRs of galaxies 
and should be viewed as complementary to dynamical mod- 
els describing the morphological transformations eventually 
related to these processes. Both starbursts and SF trunca- 
tion/strangulation require additional free parameters. The 
first of these is the time of onset tb, i.e. the age at which a 
previously undisturbed galaxy begins to feel changes in its 
SFR. Starbursts are described by a prompt increase in SFR 
followed by an exponential decline on some given timescale 
■7%, that typically is of the order of a few hundred Myr for 
normal-size galaxies and of order 10 5 — 10 6 yr for dwarf 
galaxies. 

Across the literature, a number of different definitions 
for the strength of a burst are in use. The strength of a burst 
is either described by the amount of stars formed during the 
burst as compared to the stellar mass at the beginning or at 
the end of the burst or by the amount of all gas available at tb 
that is transformed into stars during the burst, equivalent to 
the star formation efficiency (SFE) of the burst. We prefer 
the latter definition and define the burst strength via the 
star formation efficiency of the burst: 

„ stellar mass formed during burst 

SFE burst) = — (7) 

gas mass available at onset of burst 

This definition limits burst strengths to be within the 
range ... 1 and gives, at the same time, a measure of the 
global SF efficiency during the burst. Please note that since 
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the burst strength depends on the amount of gas available 
at the onset of the burst, bursts of the same strength but 
occuring at different times and/or in galaxies of different 
types do not necessarily form the same amount of stars. 
E.g. a strong burst in an old and gas-poor Sa galaxy will 
form a much smaller mass of new stars than a weak burst 
early in the life of a gas-rich Sd galaxy. All burst strengths 
can only be determined at the end of the bursts. Before the 
end of the burst only lower limits can be estimated. 

A truncation or strangulation of the SFR, e.g. as a 
galaxy falls into a galaxy cluster and has its Hf stripped by 
ram pressure, is described either by an abrupt termination 
of SF or an exponential decline of the SFR on a timescale t&. 
After a burst or after SF truncation/strangulation, GALEV 
models can have SFRs going either exponentially to zero or 
to some constant low level. 



6 THE WEB INTERFACE 

We present the new web interface that enables everyone to 
run her or his model of choice in a fast, easy to learn and easy 
to use, and comfortable way over the internet. The interface 
can be found at 

http: / /www. galev.org 

Using a web interface enables us to impose some checks 
previous to program execution dealing with the most seri- 
ous and frequently occuring problems. For the user this has 
the important advantage of always having the latest version 
at hand, including the most up-to-date input physics like 
isochrone sets and stellar libraries. 



6.1 How to use the web-interface 

The web-interface consists of four steps: Two for the input 
of the model parameters and the requested output, one for 
parameter checking and at last the actual running of the 
program on the web-server. 

6.1.1 Principal model parameters 

On the first page the user has to decide which type of model 
(s)he wants to run. Those defining parameters are 

(i) Galaxy type: choose between the most common 
types E, SO, and Sa to Sd for undisturbed galaxies, for which 
the parameters described in Sect. 15.21 will be used. In addi- 
tion we offer free types, defining the shape of the SFH as de- 
scribed by eqs.[4H6] We also offer instantaneous burst (SSP) 
models, as well as completely definable SFHs. For the latter 
mode the user has to specify a file with SFRs as function of 
time in the next step. 

(ii) Burst: we currently offer three possibilities: No burst, 
i.e. an undisturbed galaxy, a Burst with given strength, du- 
ration, and exponentially declining SFR, beginning at some 
time tb and SF Truncation with the SFR declining exponen- 
tially to zero on some specified timescale, again from time 
tb on. 

(iii) Number of filters for which magnitudes are to be 
computed. 



If the user chooses to run a standard model for an undis- 
turbed galaxy, i.e. an elliptical or spiral galaxy, then GALEV 
in the following uses the parameters described in Sect. 1531 
For all other cases, i.e. either user-defined galaxy types or 
models with burst /truncation, the parameters will be en- 
tered in the following step. 

6.1.2 Model configuration 

Depending on the choices from the previous page the user 
might not have all possible parameters to choose from; only 
the ones required for the specified model type will be shown. 
Those are: 

(i) We off er a set of Ini tial Mass Functio ns (IMF), 
namely from lSalpeterl IT955h andlKroupal (|200fl ). both using 
the Lejeune (|Leieune et aljl997l.ll99a) stellar a tmosphere li- 
brary. Further IMF shapes (e.g. lchabrierll2003f ) and different 
sets of stellar libraries are in preparation. In a future release 
we plan to also allow the user to choose customized IMFs; 
this feature, however, will be limited to SSP models (i.e. for 
star clusters). 

(ii) Gaseous emission can be switched off or on, and the 
user can choose to include only continuum emission or both 
continuum and line emission. In fixed metallicity models, the 
gaseous emission is evaluated for the respective metallicity, 
in the chemically consistent case, it is calculated appropri- 
ately as described in Sect. 13.31 

(iii) The metallicity: choose between chemically consis- 
tent models, i.e. those including all enrichment effects de- 
scribed above, or a model with metallicity fixed to one of 
the following values [Fe/H] = -1.7, -0.7, -0.4, 0.0, +0.4. 

(iv) The galaxy type is shown as a reminder, but cannot 
be changed here any more. To change the type the user has 
to go back to the first step. 

(v) Galaxy mass is another free parameter, mainly used 
for normalization. We do not include a mass-metallicity re- 
lation of any kind, so this parameter has no impact on the 
resulting colours. Note that when GALEV is combined with a 
cosmological model to e.g. obtain apparent magnitudes, the 
model galaxy mass can optionally gets rescaled to match, 
after a Hubble time, the average B-band luminosity of local 
galaxies of the respective type in Virgo. 

(vi) SFH parameters: if (and only if) one of the free 
models has been chosen, then the variables a (normalization 
for the SFR of an elliptical model), r (e- folding timescale 
of an elliptical model), j3 (factor of proportion for spiral 
models), and tpo (constant SFR of an Sd model) can be 
freely chosen by the user. The meaning of those parameters 
is explained in more detail in Sect. 15.21 If in the previous step 
the user requested the SFH to be read from a file this will 
be specified here. We caution the reader to carefully study 
the web-output from GALEV, since depending on the chosen 
parameters there might be problems, e.g. a SFR exceeding 
the amount of gas that is available. 

(vii) If the user wants to compute a galaxy featuring a 
burst or SF truncation, then (and only then) he/she has 
to specify the time of the onset of the burst, expressed in 
years after galaxy formation and the e-folding timescale for 
the decline of the SFR during the burst. While those two 
parameters are common to both truncation and burst, only 
in the latter case one also needs to supply a burst strength. 
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Following our definition of the burst-strength in Sect. 15.31 
(eq. [7}, this specifies the fraction of gas available before the 
burst that is to be converted into stars during the burst. If 
one instead wants to specify the SFR at the very onset of 
the burst, SFR(t = tb), those two numbers can be converted 
into each other by b = SFR(t = tb) x T b /M gas (t = t b ). 

(viii) GALEV currently supports two differe nt extinction 
laws t hat can be applied to the spectra, the 
1994) law for starburst galaxies, and the 



Calzetti et al 



Cardelli et al 



1989j) Milky Way extinction law, using the standard value 
of Rv = 3.1. To compute colours for different extinction 
values, the user can specify a maximum E(B-V) value and a 
step-size. The output (see below) will then include multiple 
files, one for each extinction value in the requested range. 

(ix) Cosmological parameters need to be specified to 
model the evolution of galaxies with redshift. The user can 
choose Hubble constant Ho , Om , and £Ia , while J7k is fixed to 
0. To convert ages into redshifts the user also has to select 
a formation redshift Zf oxm , so that the age of the galaxy 
is given by age(z) = tir(z) — tH(zform) with tH being the 
Hubble time at redshifts z and Zf orm . Note that the modeling 
process is terminated at a galaxy age of 16 Gyr. A choice of 
cosmological parameters leading to a galaxy age of more 
than 14 Gyr (e.g. too low a Hubble constant) hence results 
in an error-message, so that no colours will be computed for 
this case. 

(x) Although GALEV offers a full range of different output 
options (see below), not all the numbers will be actually 
needed for any specific application. The user can therefore 
use the section Output parameter to restrict the output, 
leading to faster execution and smaller downloads. 

(xi) To compute magnitudes from the spectra the user 
can choose from a large list of filter functions and also 
specify the magnitude system on which magnitudes are to 
be based (Vega, AB or ST magnitudes) For each filter a dif- 
ferent magnitude system can be chosen. This eliminates the 
uncertainties involved in any a posteriori transformations 
between different systems. 

We also offer an extensive online help giving examples 
for different parameters and how they affect the resulting 
SFH. 



6.1.3 Output 

In the following we will describe which outputs are available 
and also how they are computed. 

(i) The integrated spectra for each timestep are the 

most important output, since from those all magnitudes and 
redshifted spectra are derived. 

(ii) Convolving the spectra with the appropriate filter 
functions and applying the appropriate zeropoints gives us 
the magnitudes of the model as a function of time. 

(iii) GALEV also delivers the most useful diagnostic data 
for each timestep, such as current stellar and gaseous 
mass, star formation rates, ISM metallicity, and ion- 
izing flux NLyC- 

(iv) After applying a cosmological model to convert 
galaxy age into redshift (assuming a formation redshift) 
we can compute the redshifted galaxy spectrum, that in- 
cludes all evolution ary effec t s. We also apply the intergalac- 
tic attenuation from lMadaul (1 19951 ) to account for absorption 



of flux shortwards of Lya by intervening neutral hydrogen 
clouds. 

(v) Convolving those redshifted spectra with filter func- 
tions FF(A) and adding the bolometric distance modulus 
gives us apparent magnitudes with and/or without attenua- 
tion, that then can be used e.g. for comparison with observed 
SEDs to derive photometric redshifts or ages and metallici- 
ties. 

(vi) From the data described above we can derive cosmo- 
logical corrections due to the shifting of the filter functions 
to shorter restframe wavelengths (k-corrections) and evo- 
lutionary corrections (e-corrections). Those corrections are 
computed as follows: 



k(z) = -2.5 x log 
e(z) = —2.5 x log 



/ ao f(t = t 0j z,A) x FF(A)dA 
/ °°f(t = t 0) 0,A) x FF(A)dA 

XT f (* = t( z )' z - A ) x FF ( A ) dA 
j oo f(t = t ,z,A) xFF(A)dA 



(8) 
(9) 



where f(t, z, A) is the flux at wavelength A of a galaxy of age 
t at redshift z, to the current age of the universe (which de- 
pends on the specified cosmological parameters), and FF(A) 
a filter function. 



6.1.4 Parameter checking and execution 

As a third step in the modeling process we perform a quick 
check to ensure that all required parameters are given and 
correspond to valid combinations. The page also displays all 
given parameters to allow the user to check the input, and 
eventually perform corrections by going back to the previous 
page. However, the tests performed are just basic validity 
tests and do not ensure that the computed model makes 
physical sense. 

In the fourth and last step, after the user has made 
sure that all input is correct, we create all files necessary for 
the actual modelling. The execution of the GALEV program 
takes several minutes to run. We urge all users to carefully 
read through the given output to check if everything ran 
smoothly. 

Directly after GALEV has created all spectra, it com- 
putes the magnitude evolution as a function of time or red- 
shift in all filters requested, again taking a few minutes to 
complete. 

For compatibility reasons all output is given as human- 
readable ascii-files, with values aligned in columns that are 
separated by spaces, so that the files can be easily anal- 
ysed and/or plotted. The resulting files are then automat- 
ically combined and compressed into one archive (.tar.gz) 
file, that can be downloaded. Each archive also contains a 
small ReadMe file listing the content of each file. The mean- 
ing of individual columns are specified at the beginning of 
each file. 



6.2 Upcoming features of the web-interface 

In its current version the web-interface supports the most 
frequently used GALEV features like computing spectra and 
colours. Further features, like the computation of Lick in- 
dices or colour magnitude diagrams are currently in the pro- 
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cess of being adapted and implemented into the webpage and 
will be accessible online in the near future. 



7 FUTURE PROSPECTS 

In the near future (Weilbacher et al. 2009, in pr e p.) w e will 
offer further s tellar libraries, e.g. iMunari et all (|2005l ) and 
ICoelho et aP <|2005h . to be able to compute high-resolution 
spectra for comparison with modern spectroscopic surveys. 
Those libraries have higher spectral resolution and shorter 
wavelength co v erage but ar e not colour-corrected in the way 
iLeieune et al.l (1 19971 . Il998l ) did for the Kurucz spectra to 
reproduce the observed colours of stars from the UV through 
the NIR over the full range of effective temperatures. 

One drawback of current galaxy evolutionary synthesis 
models is that they do not include a self-consistent treat- 
ment of dust absorption and reemission. One factor con- 
tributing to this difficulty are geometric effects and depen- 
dencies of dust masses and properties on gas content, chemi- 
cal abundances and eventually even radiation field. Early at- 
tempts to consistently include dust absorption as a function 
of gas content and [Fe /H] in collaboratio n with D. Calzetti 
were encouraging (cf. iMoller et al1l200ll ) and showed that 
every undisturbed galaxy goes through a phase of maxi- 
mum extinction of E(B — V) ~ 0.4. The redshift of this 
maximum E(B — V) phase is determined by the interplay 
between decreasing gas content and increasing metallicity. 
The predicted values at high-redshift z ~ 3 (equivalent to 
a galaxy age of 2 Gyr) agree well with obs ervations from 
ISteidel et all (| 19991 ), IShaplev et al.l (|200ll ) and lColbert et al.1 
|2006l ). 

Another aspect is the coupling of GALEV to a dynamical 
model (SPH + N-body + SF + feedback) to cope with in- 
creasingly available data from Integral Field Units allowing 
spatially resolved spectroscopy. Our early attempt to couple 
GALEV models for single stellar generations with a cosmo- 
logical structure formation simulation by M. Steinmetz was 
encouraging (cf. IContardo et al.lll998l ) and showed that this 
approach is feasible. The simulated HST images of a galaxy 
at different redshifts showed fairly good agreement with ob- 
servations but ultimately failed to reproduce the correct lo- 
cal disk sizes and parameters. A better and more detailed 
description of feedback seems to be required. This can only 
be obtained from an extensive comparison between model 
results and resolved galaxy observations. The key issue is to 
have a correct criterion for SF and a correct description of 
feedback on the relevant scales and over the full range of 
SF activity - from the lowest levels in the farthest outskirts 
of galaxies to the highly clumped and clustered SF in the 
strongest starbursts. 



8 APPLICATIONS 

GALEV models have a wide range of applications from star 
clusters and resolved stellar populations of nearby galaxies 
through integrated properties of galaxies up to the highest 
redshifts. A fair number of them have been explored so far, 
many of them hand in hand with refinements or special fea- 
tures added to the models. Here we briefly recall a few of 
them to illustrate the various features of GALEV. 



8.1 Star clusters 

The simplest stellar systems to study with GALEV are star 
clusters, so-called simple stellar populations (SSPs) where 
all stars are formed essentially within one timestep and with 
the same chemical abundances. GALEV models describe the 
time evolution of SSPs with different metallicities, including 
the gaseous emission during early evolutionary stages and 
as appropriate for their respective metallicity. They can also 
incorporate extinction within the clusters' pa rent galaxy on 
the b asis of empirical extinct ion laws from ICalzetti et al.l 
1 19941 ) or lCardelli et~ai1 (|l989l ). 

Using our AnalySED tool (| Anders et al.ll2004b1 lah. we 
can compare observations in widely spaced broad-/medium- 
band filters to a grid of GALEV models and derive physical 
cluster parameters such as ages, masses, metallicity, and ex- 
tinction between the cluster and the observer If one of the 
parameters can be externally constrained (e.g. dust-free en- 
vironment, or metallicity previously determined from spec- 
troscopy) observations in at least 3 bands are required, oth- 
erwise at least 4 bands are needed. As a large number of 
star clusters can usually be covered by a single set of ob- 
servations, this is a very efficient way to study statistically 
significant cluster samples. 

As shown both in studies based on artificial star clus- 
ters and on star clusters with ages and metallicities de- 
rived independently from CMDs, accuracies in age de- 
termination of Aage/age ^ 0.3 and in metallicity deter- 
minations of A[Fe/H]/[Fe/H] ^ 0.2 are achievable, prefer- 
ably if both a short- wavelength band (U or B) and a 
NIR-band are included ([Anders et aT]|2004al ; Ide Griis et al.l 
l2003al ; Ide Griis fc Anders! 120061 ). The U-band is crucial for 
accurate age-dating, and a NIR-band (H or K) for accu- 
rate abundances. W e successfully applied these models in 
lAnders et al.l (|2004bl ) to the interpretation of the young star 
cluster systems in th e starburst dwarf galaxy NGC 1569, in 
lAnders et all (|2007l) to the analysis of the star clusters in 
the interacting Ant ennae galaxies (NGC 4038/39), and in 
iKotulla et al.l (|2008l) to the derivation of ages and metallic- 
ities of the globular clusters in the Virgo SO galaxy NGC 
4570. 



8.2 Colour-magnitude diagrams 

One special feature of GALEV is its ability to compute colour- 
magnitude diagrams (CMDs) in any desired passband com- 
bination. This is possible not only for SSPs or instantaneous 
bursts, but also for composite stellar populations with com- 
plex SFHs. This has successfully been used to identify the 
best possible passband combination to disentangle age and 
metallicity effects in star clusters in various age ranges (cf. 

iFritze et al.ll2006l). 

In lFritze fc Lilly! (|2007l ) we compared the SFH obtained 
from the CMD with those obtained from the integrated spec- 
trum, from Lick index measurements and from multi-band 
photometry in their respective accuracies and limitations. 
The basic result was that none of the methods allows to look 
back beyond a recent burst or some recent phase of enhanced 
SFR, and that all methods face very similar accuracy limi- 
tations at look-back times beyond ~ 1 Gyr. Only within the 
most rece nt Gyr, CMD analysis achieves the most detailed 
SFHs (cf. iLillv fc Fritzell2005bl lah. The importance of this 
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kind of comparative investigation lies in the fact that CMD 
analyses can only be done for the resolved stellar popula- 
tions within the nearest Local Group galaxies. All attempts 
to explore the SFHs of more distant galaxies have to rely on 
integrated spectra and, for the most distant ones, on inte- 
grated photometry only. 



8.3 Undisturbed galaxies 

In iBicker fc Fritzd (|2005l ) we used galev models to study 
the effects of the chemical evolutionary state of galaxies on 
their star formation rate indicators (Ha, [Oil], NUV and 
FUV luminosities) and found that all of them significantly 
depend on metallicity, with errors in the worst cases of up to 
facto rs of a few, confirm i ng previous observational evidence 
(e.g.. IJansen et ail 1200 ll ; iHopkins et"ai1l2003l ; iKewlev et ail 

l2004h . 

In lSchufJ et al.l (|2003h study the time and redshift evo- 
lution of bulge-to-disk light ratios in different wavelength 
bands by assuming a short timescale for SF for the bulge 
component *l/buigc ~ exp(— 1/1 Gyr) and a constant SFR for 
the disk component ^disk ~ constant. The integrated spec- 
tral and photometric evolution of different spiral galaxy 
types was then obtained by adding up the bulge and disk 
components in mass ratios so as to give, after a Hubble time 
of evolution, the observed average B-band bulge-to-disk light 
ratios for the respective spiral types Sa through Sd. This 
study showed a significant wavelength dependence of the 
bulge-to-disk light ra tios in agreement with observations by 
lEskridge et al.l (2002J). This has implications for galaxy clas- 
sification in different redshift intervals. It also opens a new 
possibility to explore bulge formation scenarios and bulge 
formation redshifts by comparing bulge-to-disk light ratios 
measured in different bands. 



8.4 Interacting galaxies and mergers 



8.3.1 Starbursts in Blue Compact Dwarf Galaxies 

Starbursts were first investigated with galev models in the 
context of Blue Compact D warf Galaxies (BCDGs) in a se- 
ries of papers by H. Kroger (IKriiger et al.lll99ll . Il992l . 1 19931 ; 
iKriiger fc Fritzel 19941 ; IKriiger et al.lll995T l. These authors in- 
vestigated a sample of BCDGs with optical and NIR pho- 
tometry in order to derive their burst strengths and the age 
of their underlying stellar population. The main results we 
want to recall here are that even very weak ongoing bursts 
can completely dominate the light in the optical, in par- 
ticular at the low metallicities Zq/50...Zq/5 typical for 
BCDGs. An underlying old galaxy component can only be 
detected in the NIR and was found for every BCDG of 
our sample. Very accurate age-dating was possible for those 
BCDGs which showed a 4600 A bump caused by WR-stars in 
their spectra. Burst strengths were found to be of the order 
of a few percent only, when defined in terms of stellar mass 
increase. They were also shown to systematically decrease 
with increasing galaxy mass, where the latter included the 
important mass contributions of HI. This result is in agree- 
ment with expectations on the basis of t he stochastic self- 
propa gatin g SF scenario put forw ard by iGerola fc Seider] 
(| 19781 ) and lSeiden fc Gerolal (|l979h . 



In lFritze fc Gerhard! (|l994al lbh we studied a grid of starburst 
models with bursts of various strengths occurring in Sa, . . ., 
Sd spirals at different ages in their spectral, photometric and 
chemical evolution and then analyzed the starburst in the 
gas-rich massive spiral-spiral merger NGC 7252. We found 
this burst to have started about 600 — 900 Myr ago and to 
have been stronger by 1 — 2 orders of magnitude than those 
in BCDGs. The bulk of information available for this galaxy 
even allowed us to estimate the SF efficiency on the basis of 
a comparison of the stellar mass formed during the burst. 

Our conservative estimate for the overall SF efficiency 
(see equation [7} during this interaction- induced starburst 
indica ted a very high value SFE ^ 0.35 |Fritze fc Gerhard] 
Il994al lbl; iFritze fc Burkertl Il995h . again about two orders 
of magnitude higher than any SFE measured for molecu- 
lar clouds in the Milky Way or the Magellanic Clouds, and 
high enough to allow for the formation of a new generatio n 
of globular clusters (cf. iBrown et all 1 19951 ; iLi et al.ll2004r i. 
GALEV models indicated that NGC 7252 at present still fea- 
tures a low-level ongoing SFR of ~ 3Mq yr _1 in its centre, 
powered for ~ 50 % by gas set free at present from dying 
burst stars and for ~ 50 % by HI falling back onto the main 
body from the tidal tails. Both of these gas delivery rates 
will decrease over the next 1 — 3 Gyr. Depending on whether 
the SFR will cease completely or continue at some very low 
level, NGC 7252 will spectrally evolve into an elliptical or 
SO galaxy over the next 1 — 3 Gyr. Already at present, NGC 
7252 features an r 1//4 — light profil e across a radial range of 
~ 14 kpc, if azimuthally averag ed (|Schweizerlll98l ) 

Independent confirmation for the unexpectedly high 
value found for the SFE in NGC 7252 came from the de- 
tection of a rich population of massive compact star clusters 
with ages in agreement with the global starburst age which, 
in turn, is in agreement with dyn amical merger ages from 
N-body + SPH simulations (e.g. Iffibbard fc Mihosl Il995h . 
Spectroscopy of the brightest clusters confirmed their metal- 
licities to be between Zq/2 and Zq, as expected if they 
formed out of the gas pre-enriched in Sc-type spirals - with 
some evidence for a moderate amount of self-enrichment 
during the burs t on the basis o f their slightly enhanced 
[a/Fe] ratios. In lFritze fc Burkertl (|l995T ) we estimated that 
the number of clusters with masses in the range of Galactic 
globular cluster (GC) masses that formed in the burst and 
survived until the present is of the same order of magnitude 
as the number of GCs present in two average-luminosity 
Sc-type spirals before the merger. Hence, this spiral-spiral 
merger will, after the fading of the post-starburst and af- 
ter the fading and dissolution of the tidal features, evolve 
into an elliptical or SO galaxy with a normal GC specific fre- 
quency. GCs of age 0.5 ... 1 Gyr have already survived the 
most critical phase in their lives, the infant mortality and 
early mass loss sta ges, and stand fair chances to survive for 



many more Gyr (cf. Lamers et alj|2005l ; iBastian fc Goodwin! 



_mv 

120061 ; iParmentier fc Fritzell200gh 

In lde Griis et al.l (|2003bh we used GALEV SSP models to 
analyze the luminosity-weighted ages pixel-by-pixel on ACS 
images of the interacting galaxies Tadpole and Mice and 
studied their star cluster populations. A surprising result 
was that about 35 % by mass of all recent SF went into the 
formation of star clusters in both galaxies, not only across 
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the main bodies of both galaxy systems but all along their 

very e xtended tidal tails. 

In iTemporm" fc Fritzel I 20061 ) we applied galev mod- 
els to investigate the SF and starburst histories of galaxies 
in a very compact group of galaxies on the basis of multi- 
band photometry and spectra and in IWehner et al.1 (|2006l ) 
we studied the SF activity and its history in the extended 
tidal debris surrounding the starburst galaxy NGC 3310. 



8.4-1 Tidal Dwarf Galaxies 

Not only star clusters can form in the low-density environ- 
ments of tidal tails but sometimes even star-forming ob- 
jects with masses in the range of dwarf galaxies: so-called 
Ti dal Dwarf Galaxies (TDGs ), or better TDG cand i dates . 
In Weilbacher fc Fritz 1 (|200ll ) and IWeilbacher et al.l (|2002l. 



l2003bl 



we analyzed the first reasonably sized sample of 
TDGs and found that they all contain a stellar population 
inherited from the spiral disk out of which the tidal tail has 
been torn, together with a significant young stellar popula- 
tion that must have been formed in situ within the tidal tail 
after it had been ejected. A characteristic feature of TDGs 
is that they do not follow the luminosity-metallicity relation 
of dwarf galaxies but all have similar metallicities character- 
istic of the HII region abundances in spiral disks. 

Again it turned out that optical observations alone are 
not sufficient to disentangle the mass contributions of the in- 
herited versus the starburst components. Even a 90 % mass 
fraction in the inherited component can be entirely hidden 
in the optical by an ongoing burst that only makes up for 10 
% of the mass. Only in optical-NIR colours can the inherited 
component be detected that is not entirely old but contains 
the mix of stellar ages present in the disk before the tidal 
tail was thrown out. 



8.5 Galaxy transformation in groups and clusters 

A variety of scenarios are discussed in the literature to ex- 
plain the transformation of the spiral-rich field galaxy pop- 
ulation into the S0-/dSph-/dE-rich galaxy population ob- 
served in rich galaxy clusters at low redshift. The follow- 
ing processses have been proposed: High-speed disruptive 
galaxy-galaxy interactions called harassment, interactions 
between galaxies and the dense hot Intracluster Medium 
(ICM), and enhanced merging within infalling groups. All 
of these scenarios are observed to be at work in a number of 
individual cases. Their relative importance, their timescales, 
transition stages, and end products, however, are not known 
yet. All of these transformation processes both affect the 
morphological appearance of galaxies and - via their SF 
histories - their spectral properties. How the timescales 
for morphological transformation and spectral transforma- 
tion relate to each other in the various scenarios and en- 
vironments is not clear to date. Removal of gaseous halos, 
outer, and inner HI disks leads to SF strangulation on long 
timescales or to SF truncation on shorter ones. Destabi- 
lization of disks through encounters or shocks as well as 
mergers within infalling groups may lead to starbursts. We 
explored aspects related to the spectral transformation of 
galaxies through all these scenarios and investigated which 
scenarios in which type of progenitor galaxy and at which 



evolutionary stage can lead to the observed luminosity and 
colour ranges of SO galaxies by implementing SF stran- 
gulation/truncation on different timescales wi th and with- 
out pr ece ding starbursts into gale v models in lBicker et al.l 
|2002l ). In lFalkenberg et all (|2009al lbl) we extended the mod- 
els to also include the evolution of the D4000 and H<5 Lick 
indices into the GALEV models for galaxy transformation and 
investigated under which conditions t he so-called E+A-, or 
k+a-, and the H^-s trong galaxies (cf. iPoggianti et al 1 12004 
iDressler et~aT]|2004 ) are formed, what is the lifetime of this 
respective phase, what is the colour and luminosity of the 
galaxy in this transition stage and what is the end-product. 



8.6 High redshift galaxies and photometric 
redshifts 

If coupled to a cosmological model, GALEV can be used to 
study the evolution of galaxies from the very onset of SF 
in the early universe until today. Accounting for the signifi- 
cantly sub-solar metallicities observed at high redshifts - in 
particular when dealing with intrinsically faint galaxies that 
dominate in deep field surveys - allows us to determine more 
accurate photometric redshifts, as comp ared to those ob- 
taine d with solar-metallicity models only (|Kotulla fc Fritzel 
120091 ) or observed temp l ates ( Kotulla 2009, in preparation). 

In iFritze fc Bicker! (|2006l ) we examined starbursts and 
their respective post-starburst stages across a wide range of 
redshifts with the surprising result that dust-free models in 
their long post-burst phases after strong starbursts at high 
redshifts can get the colours and luminosities of Extremely 
Red Objects (ERQs) as e.g. observed in the K20-survey (cf. 
iDaddi et aLll2002l ; ICimatti et aill2002f ). 

Only when we also include - in addition to our set of 
undisturbed models E, Sa, . . ., Sd - the very blue starburst 
phases and also their extremely red postburst phases can we 
reproduce the full range of colours observed e.g. in the Hub- 
ble Deep Fields (Kotulla & Fritze 2009, in preparation). In 
Fig. [9] we show the F606W-Ks (approx. V-K) colour evo- 
lution for a set of undisturbed models E, and Sa through 
Sd, and for models with major starbursts (burst strengths 
chosen to consume 70% of the in each case available gas) 
occurring at different ages in a previously undisturbed Sb 
galaxy. A comparison to th e photometric galaxy catalog of 
iFernandez-Soto et all (| 19991 ) for the HDF shows that galev 
models can describe the full colour range, even without any 
dust. Note that dust certainly plays a role in ongoing star- 
bursts but not any more during postburst stages. Appar- 
ently, as observed in the local Universe, galaxies use up and 
destroy part of their dust in starbursts, while the rest is 
blown away by the end of the burst. 



8.7 Redshift evolution of ISM abundances: spiral 
models vs. DLAs 

The SFHs of our closed-box model galaxies were constrained 
by requiring agreement after 12—13 Gyr of evolution not 
only with observed spectrophotometric properties of the var- 
ious galaxy types but also with their gas content and ISM 
abundances (cf. Sect. [5j5) . While this ensures agreement at 
the present stage it needs not necessarily imply agreement 
over the entire evolutionary path. Some subtle interplay be- 
tween SFR and infall rate, e.g., could lead to the present-day 
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Figure 9. V-K colour as function of redshift for five undisturbed 
spectral types E and Sa through Sd. Extremely blue colours are 
found during active star bursts, while extremely red colours can 
be reached during the postburst phases of early starbursts, here 
shown for galaxy ages of 1, 2, 3, 5, and 10 Gyr. 



agreement after ages of disagreement. Infall rates are very 
hard to constrain from a comparison with spectrophotomet- 
ric properties alone. They can, however, be constrained by 
comparing chemical abundances. It is therefore of prime in- 
terest to compare the redshift evolution of ISM abundances 
to observations of high-redshift galaxies. 

Transforming the time evolution of ISM abundances 
into a redshift evolution only requires the one-to-one trans- 
formation between galaxy age and redshift that is given by 
any set of cosmological parameters (see. Sect. 13. 9[) 

Damped Lyman-a absorbers (DLAs) are a particular 
class of absorption systems, usually observed in the line- 
of-sight of quasars, with damping wings on the Lyman-a 
line due to the high column densities of neutral hydrogen of 
N(HI) ^ 10 21 cm~ 3 in the absorbers. The damped Lyman-a 
line is always accompanied by a number of low-ionisation 
heavy element lines, many of which are the dominant ioni- 
sation states in those absorbers. DLAs are most easily ob- 
served in the redshift range z ~ 2 to z > 4 and have been pro- 
posed to be potential progenitors of present-day spiral disks: 
Their HI column densities are similar to those in local gas- 
rich disks and their gas content, derived from rotational ve- 
locities of the order 100 . . . 200fcm s -1 , is similar to gaseous 
plus stellar masses in local spirals. In lLindner et al.l(|l999l ) 
we compared the redshift evolution of our spiral galaxy mod- 
els, calculated in a chemically consistent way with stellar 
yields for a large number of elements and 5 different stellar 
metallicities to the first reasonably sized set of Keck HIRES 
abundances for DLAs from the literature. We found surpris- 
ingly good agreement for all the different abundances over 
the entire observed redshift range with our closed-box mod- 
els. 

Studying the impact of various amounts of infall showed 
that it was not possible to accommodate more than a moder- 
ate amount of infall, increasing the total mass from redshift 
2 until the present by not more than a factor ~ 2, with- 
out losing agreement with either the spectral or the chem- 



ical properties or both, no matter how we tuned the SFH. 
We hence concluded that from a chemical evolution point of 
view DLAs might well be the progenitors of present-day spi- 
ral galaxies of all types Sa . . . Sd and that those are reason- 
ably well described by closed-box or moderate-infall models. 
This implies that DLAs have already almost all the mass of 
present-day spirals - albeit almost completely in the form of 
HI gas. GALEV models indicated low stellar masses and low 
luminosities for DLA galaxies, in agreement with the large 
number of non-detections/upper limits. Mass estimates from 
rotation velocities derived from detailed modeling of asym- 
metries detecte d in some DLA p rofiles confirmed our mass 
predictions (cf. I Wolfe et~alll2005l) . GALEV models also pre- 
dicted a change in the DLA galaxy population from high to 
low redshift and showed that the DLA phenomenon can be 
understood as a normal transition stage in the life of every 
spiral. During their enrichment process to higher metallici- 
ties they convert their gas reservoir into stars, therefore get 
increasingly gas-poor so that above a certain metallicity they 
drop out of DLA sa mples due to too low a gas content (cf. 
iLindner et"afll 19991 . for more information). This is in agree- 
ment with observations that show the lowest redshift DLA 
galaxies to be low-luminosity late-type or irregular galaxies. 



9 SUMMARY 

This paper presents the GALEV evolutionary synthesis mod- 
els for star clusters, undisturbed galaxies and galaxies with 
starbursts or/and star formation truncation now available 
on the web at 

http:/ /www. galev.org 

We describe the input physics currently used, that will con- 
tinuously be updated. 

For a number of different stellar IMFs, the spec- 
tral evolution of star clusters of metallicities in the range 
— 1.7 ^ [Fe/H] ^ +0.4 can be calculated, and a large number 
of filter systems are available for the photometric evolution 
as well as the full set of Lick absorption indices. 

GALEV features a unique combination of characteristics 
that allow for what we call a chemically consistent modelling 
of the chemical evolution of the ISM together with the spec- 
tral evolution of the stellar component. 

This means that the initial abundances of every stellar 
generation are accounted for by using input physics (stellar 
evolutionary tracks, stellar model atmospheres, gaseous line 
and continuum emission, stellar lifetimes, yields and rem- 
nant masses) appropriate for the increasing initial abun- 
dances present at the formation time of successive stellar 
generations. This chemically consistent modeling accounts 
for the observed broad stellar metallicity distributions in 
local galaxies as well as for the increasing importance of 
subsolar abundances in local late-type and low luminosity 
galaxies and in high redshift galaxies. 

Galaxy models can be calculated either in the chemi- 
cally consistent way or for some fixed metallicity upon re- 
quest. Models give spectra, emission line strengths and Lick 
absorption features, photometric quantities for a large num- 
ber of filter systems, and chemical abundances, gaseous and 
stellar masses, star formation rates, etc. in their time evo- 
lution for normal galaxies, galaxies with starbursts or/and 
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star formation truncation as specified by the user or for user 
customized star formation histories. 

If a cosmological model is selected, all quantities are 
also provided in their redshift evolution, fully accounting 
for evolutionary and cosmological corrections and including 
the attenuation by intergalactic neutral hydrogen. 

We present the models, the input physics they use, their 
calibrations, the web interface and some examples of selected 
applications for illustration and we discuss current limita- 
tions and future prospects. 
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APPENDIX A: WORKING IN PICTURES 

In the following we will present the workings of GALEV in 
a more vivid way by guiding the reader through a series of 
steps. We give examples for the input physics GALEV uses 
and show how this is used throughout the process of mod- 
elling a galaxy. 



Al General steps 

In the first step we present how one gets from isochrones, 
a stellar IMF, spectral library, and atomic physics to an 
integrated isochrone spectrum. For this example we choose 
to work with isochrones from the Padova group, a Salpeter- 
IMF from O.1M to IOOMq and stellar spectra from the 
Lejeune library, all using solar metallicity. 



Step 1: Choose set of isochrones 

Fig. IA1I shows solar-metallicity isochrones for three dif- 
ferent ages of 4 Myr, 100 Myr and 1 Gyr as a colour- 
magnitude diagram. For very young ages the main sequence 
reaches up to high masses and hence very high luminosities 
~ —10 mag. For the later stages the red giant branch (RGB) 
at (B — V) ~ 1 ... 1.4 is clearly visible and also the asymp- 
totic giant branch (AGB). Note that stars on the AGB can 
reach extremely red colours of (B — V) ~ 4 mag and at the 
same time get very bright in the NIR (Mr ~ — 10 mag) dur- 
ing the thermal pulsation (TP) phase. 
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Figure Al. Solar metallicity isochrones from the Padova group 
for three different ages of 4 Myr, 100 Myr and 1 Gyr 
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Figure A2. T wo examples of init ial mass functions, here from 
iKroupal |200ll ) and lSalpeteil [|l955l l 

Step 2: Select Initial Mass Function 

In a second step we have to choose a parameterization for the 
number of stars as function of their mass. Fig. lA2l shows two 
common IMFs, the Salpeter and Kroupa IMFs. For masses ^ 
IMq both predict roughly comparable number of stars, but 
they differ in the low mass regime ^ 1M© . This will not only 
affect mass-to-light ratios, but also have an impact on the 
resulting spectra and in particular the chemical enrichment, 
since fewer low-mass stars also means that less mass gets 
locked up in long-lived stars. 

Step 3: Populate isochrones with stars 

With both isochrones and IMF at hand we can now populate 
the isochrones with stars taken from an IMF. For each mass 
point on the isochrone the IMF tells us how many stars were 
being formed with this mass. The results are shown in Fig. 
IA3I where we also added some scatter around the individ- 
ual points to give the artificial colour-magnitude diagram 
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Figure A3. Synthetic colour magnitude diagram where each 
point from the isochrone is assigned a numbers of stars based on 
the IMF. We added some scatter around each point to produce a 
smoother appearance. 

a smoother and more realistic appearance. We also added 
at least one star to each point to bring out all phases of 
stellar evolution. However, in reality, this might not be the 
case, since some phases, e.g for very young and extremely 
hot white dwarfs last for a very short timescale only. 

Step 4: Assign a spectrum to each star 

Next we assign a spectrum to each star in the colour- 
magnitude diagram. The parameters required for this are 
either directly given in the isochrones (e.g. the effective tem- 
perature T e ff ) or can be derived from other parameters (e.g. 
the surface gravity logg can be calculated from the stellar 
mass, its bolometric luminosity and effective temperature). 
In Fig. IA4l we show four sample spectra for stars with iden- 
tical log g = 4.0, but temperatures ranging from 20000 K 
down to 2000 K. Hot stars have a smooth slope and are 
bright in the UV. Stars with T off « 10000 K have promi- 
nent Balmer absorption lines characteristic for spectral type 
A stars. Cool stars are dominated by broad molecular ab- 
sorption bands and only emit at long wavelengths. Not all 
stellar parameters are directly covered by the library, all 
other values have to be interpolated between neighbouring 
points in the library. Is is therefore crucial for the stellar 
library to cover the full parameter range. 

Step 5: Integrate spectra of all stars 

To create stellar isochrone spectra we now integrate the light 
of all stars. Special care has to be taken to account for the 
different luminosities of stars at their respective stages of 
evolution. In Fig. IA5l we show spectra for the same ages dis- 
cussed above, 4 Myr, 100 Myr, and 1 Gyr. The youngest 
spectrum is almost completely dominated by the hottest 
and brightest high-mass stars that also dominate part of the 
100 Myr spectrum. For this spectrum the highest-mass stars 
have already exploded as supernovae and with them the flux 
in the FUV already has decreased considerably. In the op- 
tical the 4000A break starts to appear, becoming stronger 
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Figure A4. Sample spectra from the Lejeune stellar library for 
stars with logg = 4.0 and T cff = (20000, 10000, 5000, 2000)K. 
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Figure A5. Integrated starlight for the three isochrones shown 
in Fig. EH 



with time. After 1 Gyr there is little UV flux remaining, 
the spectrum is now dominated by stars with lower masses 
w 2M Q . 



Step 6: Compute gaseous line and continuum emission 

Until now we have only dealt with stellar light. However, as 
mentioned earlier, at young ages isochrones contain large 
contributions from high mass stars. These emit a large 
fraction of their light in the UV and also produce signifi- 
cant numbers of photons energetic enough to ionize hydro- 
gen. The energy injected into the surrounding interstellar 
medium produces both emission lines and also continuum 
emission. The detailed mechanisms of this are discussed in 
Sect. 13.31 and shown in Fig. IA6I Note that while the strength 
of the heavy element emission lines depends on the metal- 
licity of the gas, the strength of the continuum emission and 
hydrogen emission lines purely depends on the ionizing pho- 
ton flux. 
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Figure A6. Gaseous emission spectrum for solar metallicity. 
Shown in blue is the continuum emission, and continuum and 
line emission in red. 
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Figure A7. Integrated isochrone spectra for the same isochroncs 
shown in Figs. IA1I and IA5I including both stellar and gaseous 
emission. 



Step 7: Add gas emission to stellar spectrum 

For each isochrone we now compute the gaseous emission 
spectrum based on the metallicity and ionizing photon flux 
derived from isochrones and IMF. Since the latter heavily 
depends on the age of the isochrone so does the strength of 
the gaseous emission. In Fig. IA7I we show the same spec- 
tra as in Fig IA5I but now with gaseous emission included. 
For the first spectrum at an age of 4 Myr the emission lines 
clearly stand out in the optical. However, for the later stages 
there are no changes compared to the purely stellar spectra, 
since only stars with masses > 20Mq and hence very short 
lifetimes < 10 Myr produce the majority of the ionizing flux. 

In GALEV those previous steps are repeated for all 
metallicities and ages available from the isochrone set, yield- 
ing a full set of isochrone spectra including both stellar and 
gaseous emission. 
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Figure A8. Example for contributing isochrone spectra during 
interpolation in age. 

A2 Additional steps for star clusters 

Continuing from the library of isochrone spectra GALEV can 
now compute the finer grid required for the study of star 
clusters. 

Step 8: Interpolate between ages given by the isochrones 

The time resolution offered by the isochrones in general only 
spans a relatively coarse grid with logarithmic age spacing. 
However, it is advantageous to create a grid with smaller and 
linear age steps. This can be done by interpolating additional 
ages to fill the gaps in the isochrone grid. Since spectra vary 
roughly linearly with logarithmic time, i.e. changes are larger 
at small ages and small at large ages we use this as basis for 
our algorithm. In Fig. IA8I we show two isochrone spectra 
for ages of log(t) = 6.7 and log(t) = 6.9 and the result- 
ing spectrum for log(t) = 6.85. The resulting spectrum was 
computed by Spec(log(t) = 6.85) = 0.25 x (Spec(log(t) = 
6.7) +3 x Spec(log(t) = 6.9). 

To obtain a full grid to be used e.g. for age-dating star 
clusters one has to repeat this procedure for each required 
time step and each metallicity. The steps involved with ap- 
plying dust extinction and computation of magnitudes from 
spectra is explained in Steps 11 and 14 below. 

A3 Additional steps for galaxies 

Galaxies are completely different from star clusters since 
they contain multiple stellar populations, generally covering 
a range of ages and metallicities. We explain the basic steps 
dealing with multiple populations based on a toy-model of 
a galaxy made from only two populations with different 
ages and metallicities. The subsequent explanations use the 
model of an undisturbed elliptical galaxy. 

Step 8: Compute initial abundances for next generation 

To derive the chemical evolution of galaxies we have to know 
four things: 1) The total mass of the galaxy, including both 
the mass of stars and gas. 2) The star formation history 
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Figure A9. For our toy-model we construct a very primitive 
star formation history of 2 bursts, with SFRs of 100 M© yr and 
50 Mgyr 1 , each lasting 10 6 years. 



(SFH), i.e. how many stars are formed at each particular 
time. 3) The life times of the stars formed. This in combi- 
nation with the SFH yields some form of star death record. 
And we have to know 4) the end products of each star, i.e. 
the mass of its remnant, and the mass of gas and heavy 
elements returned to the surrounding ISM. For each time 
we then have to keep track of the masses of stars, gas and 
metals in the gas. All those quantities are changed by star 
formation and the return of gas and heavy elements from 
dying stars. The ratio of metals to gas is the crucial fac- 
tor since it determines both the spectra and lifetimes of the 
newly formed stars. 



Step 9: Interpolate between ages AND metallieities 

For all galaxies with an extended SFH we will have stel- 
lar populations not agreeing with the coarse grid given by 
the isochrones. We therefore have to interpolate between 
isochrone spectra of different ages and also different metal- 
lieities. The details of this process are described in the con- 
text of star clusters above and are also shown in Fig. IA8I 



Step 10: Add up spectra weighted with SFH 

To ease the understanding how GALEV assembles a galaxy 
spectrum from the individual spectra of each of its con- 
stituent populations, we first consider a toy model of a 
galaxy made from only two populations. Both populations 
are described by intervals of 1 Myr duration each, occuring 
at an age of 100 and 200 Myr and forming stars at a rate 
of lOOMoyr" 1 and 50M©yr _1 , respectively (see Fig. Eg} . 
Both intervals are short compared to the age of the youngest 
isochrone so that they each can be described as a population 
of one age. We further assume that the earlier population 
(Burst 1) is formed with a metallicity of [Fe/H] = —1.7 
or 1/50Z©, and the second (Burst 2) with a metallicity of 
[Fe/H] = -0.7 or 1/5Z©. 

We will now show how to derive the spectrum of our toy- 
galaxy at two ages of 300 Myr and 1 Gyr. For each timestep 
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Figure A10. Interpolated isochrone spectra for each burst and 
the resulting galaxy spectrum, for galaxy ages of 300 Myr (top 
panel) and 1 Gyr (bottom panel). 



t one needs to derive the ages r and metallieities of all pop- 
ulations formed prior to this time. Those are then weighted 
by the SFR at their respective formation time SFR(t — r) 
multiplied with the length of a time step. Those are then 
added up to yield the galaxy spectrum. For our toy galaxy 
at an age of 300 Myr we know that Burst 1 has metallicity 
[Fe/H] = -1.7 and an age of 200 Myr. Burst 2 has metal- 
licity [Fe/H] = -0.7 and an age of 100 Myr. The resulting 
spectrum thus can be constructed by added up those two 
isochrone spectra. The resulting spectrum is shown in the 
upper panel of Fig. IaToI 

The spectra for all other times are created in a similar 
way. Metallieities for each population stay the same, while 
ages increase with time. The resulting spectrum of our toy- 
model galaxy at an age of 1 Gyr is shown in the bottom 
panel of Fig. IaToI 

In the following we will leave our toy-model galaxy and 
instead show the remaining steps that are necessary to derive 
magnitudes for an elliptical galaxy with a small amount of 
dust at redshift z = 3. 
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Figure All. Relation between galaxy age and redshift for differ- 
ent cosmological parameters and in particular varying formation 
rcdshifts. 
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Figure A12. Extinction curves from Calzetti and Cardelli. Both 
show very little transmission at short wavelengths indicating that 
most of the light gets absorbed. 



Step 11: Apply evolutionary correction 

Since we want to model the galaxy at a cosmologically sig- 
nificant redshift of z = 3 we have to take evolutionary cor- 
rections into account, i.e. we observe the galaxy at an earlier 
evolutionary state. We therefore have to know the age of the 
galaxy at this redshift. In Fig. IA11I we show the redshift- 
galaxy age relation for a small range of cosmological param- 
eters. It is obvious that the density parameters Qm and Qa 
influence the solution, but also the formation redshift Zf orm , 
the redshift at which the galaxy started forming stars. The 
impact of those evolutionary corrections can be seen in Fig. 
IA14I The first row shows the galaxy at redshift z=0 with 
an age of ~ 13 Gyr. The second row shows the galaxy at its 
evolutionary state at z — 3; the galaxy at redshift 3 started 
forming stars only sa 1.5 Gyr earlier. 

Step 12: Apply extinction 

In a next step we apply the attenuati on due to inte r stellar 
dust. For our example we choose the ICalzetti et all (|l994f ) 
extinction law (see Fig. IA12I and choose an intermediate 
degree of extinction, E(B-V)=0.2. F or comparison we als o 
show the dust attenuation curve of ICardelli etafl (|l989h . 
Both extinction curves have in common that the transmis- 
sion, i.e. the fraction of light that remains unabsorbed, drops 
towards shorter wavelengths, leading to a reddening of the 
galaxy light. The results on the spectrum can be seen from 
the difference between the second and third row in Fig. IA14I 

Step 13: Redshift spectrum 

The spectrum now has to redshifted by a factor (1+z). Note 
that also the flux has to reduced by the same factor (1+z) 
to properly account for the cosmic expansion. 

Step 14: Apply intergalactic attenuation curve 

Due to the high redshift at which we observe our galaxy we 
have to correct for intergalactic absorption due to interven- 
ing neutral hydrogen clouds. These absorb light shortwards 
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Figure A13. Transmission of the intergalactic medium as func- 
tion of observed frame wavelength for sources at different red- 
shifts. 

of the Lyman-a line (1216 A) and hence reduce the flux in 
this region. Fig. IA13I shows the transmission of the IGM as 
function of observed frame wavelength for sources at differ- 
ent redshifts z — 1 ... 6. For a galaxy at redshift 3 this means 
that w 30% of the light between the Lyman-break (912 A) 
and the Lyman-a line is absorbed, while shortwards of the 
Lyman-break essential all flux is absorbed (hence the name 
Lyman-break) . Fig. IA14I shows in the two bottom rows the 
spectrum with (5th row) and without (4th row) this atten- 
uation. 

All the previously mentioned effects are summarized in 
Fig. IA14I The first row shows the spectrum of the galaxy at 
redshift z = at an age of 13 Gyr. The second row still is 
at z = 0, but at an evolutionary state already corresponding 
to 2 = 3. The following panel shows the spectrum with a 
reddening of E(B — V) = 0.2 mag; Here most of the far-UV 
flux is already absorbed by dust. The fourth row shows the 
spectrum redshifted to z = 3 and the last row also takes 
intergalactic attenuation into account. 
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Figure A14. Galaxy spectrum of an elliptical galaxy at different 
stages of the modelling process. Each following spectrum includes 
all effects shown above. 



Step 15: Convolve with filter response curve to compute 
magnitudes 

We can now convolve the final spectrum with filter response 
curves. Therefore each point in the spectrum is weighted 
with the relative filter response at the corresponding wave- 
lengths and then integrated over all wavelengths. To convert 
the resulting fluxes into observable magnitudes we have to 
apply zeropoints according to the requested magnitude sys- 
tem, e.g. Vega or AB. In Fig. IA15l we show a wide selection 
of currently available filters from different space and ground- 
based telescopes. 
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Figure A15. Selection of the filters currently offered by GALEV, 
ranging from the far-UV to near-infrared. 
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Step 16: Add distance modulus 

In a very last step we can convert the absolute magnitudes 
obtained with the filters into apparent magnitudes. This is 
done by simply adding the bolometric distance modulus for 
the particular redshift. In Fig. IA16I we show the evolution 
of the distance modulus with redshift for a small selection 
of possible cosmological parameter sets. For the previously 
studied galaxy at redshift z = 3 we have to add a distance 
modulus of 47 mag. 

This paper has been typeset from a Tp^/ I^TgX file prepared 
by the author. 



Figure A16. Bolometric distance modulus as function of redshift 
for a range of cosmological parameter sets. 
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